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Abstract: We propose a formula for computing the (moduli-dependent) contribution of 
multi-centered solutions to the total BPS index in terms of the (moduli-independent) indices 
associated to single-centered solutions. The main tool in our analysis is the computation 
of the refined index Tr {—yY^^ of configurational degrees of freedom of multi-centered BPS 
black hole solutions in A/" = 2 supergravity by localization methods. When the charges 
carried by the centers do not allow for scaling solutions (i.e. solutions where a subset 
of the centers can come arbitrarily close to each other), the phase space of classical BPS 
solutions is compact and the refined index localizes to a finite set of isolated fixed points 
under rotations, corresponding to coUinear solutions. When the charges allow for scaling 
solutions, the phase space is non-compact but appears to admit a compactification with finite 
volume and additional non-isolated fixed points. We give a prescription for determining the 
contributions of these fixed submanifolds by means of a 'minimal modification hypothesis', 
which we prove in the special case of dipole halo configurations. 
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1. Introduction and summary 

In A/" = 2 supersymmetric string vacua, BPS states in suitable large charge limits can be 
represented as multi-centered black hole solutions oi N = 2, D = A supergravity [1, 2, 3]. To 
compute the moduli-dependent index fi(7;t) associated with such configurations, one needs 
to combine two independent sets of data. 

The first part of the data are the indices ^^'^(7) associated with single centered BPS black 
holes carrying electromagnetic charges 7. In the supergravity approximation, the index is 
given by the exponential of the Bekenstein-Hawking entropy [4, 5], and is independent of the 
asymptotic values of the scalar fields t (within a given basin of attraction) by virtue of the 
attractor phenomenon [6, 7, 8]. Effects of classical higher derivative corrections to the low 
energy effective action can be incorporated by using Wald's modification of the Bekenstein- 
Hawking formula [9, 10, 11, 12], while quantum corrections to the index can in principle be 
computed using the quantum entropy function formalism [13, 14, 15]. 

The second part of the data is the index of the supersymmetric quantum mechanics de- 
scribing multi-centered black hole configurations. In this description, the centers are treated 
as pointlike, entropy-less particles carrying (in general, mutually non-local) electromagnetic 
charges ai, . . . , a„, and kept in equilibrium by balance of forces [1]. The space of solutions of 
this mechanical problem is a 2{n — l)-dimensional symplectic space A^„ with an hamiltonian 
action of the rotation group 5*0(3), which can be quantized by the standard procedure of 
geometric quantization [2, 16]. Unlike the first part of the data, the index of these con- 
figurational degrees of freedom, which we denote by g{{ai\]t)^ depends sensitively on the 
asymptotic values of the scalar fields t. In our previous work [17], reviewed in [18], we showed 
how to compute the jump of this index across a wall of marginal stability by localization 
with respect to a U{1) subgroup of 5*0(3) corresponding to rotations along the z-axis. In 
this work, we extend the techniques of [17] to compute the configurational index away from 
the walls of marginal stability, and propose a formula to combine this result with the indices 
associated to single-centered black holes in order to compute the total BPS index. 

The approach of [17] was based on several simplifying facts. First it was shown that 
since identical centers do not interact, one can replace the Bose- Fermi statistics carried by 
the centers by Boltzmann statistics, provided the index VL{ai) carried by the centers - or 
more generally the 'refined index'^ Vt^ciioHiU) = Tr '(— y)^"'^, where J3 denotes the generator 

^Here Tr' denotes the trace over states carrying a fixed set of charges after removing the contribution 
from the fermion zero modes associated with broken supersymmctry generators. While Tr'(— y)^''^ is in 
general not a protected index away from y = 1 (except in rigid Af ^ 2 field theories), it is essential to 
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of rotations along the z axis and ?/ is a real parameter - is replaced by an effective rational 
index firef (cti, y) [19, 20, 17, 2f ]. Thus, even when some of the a^'s are identical, one may still 
treat the centers as distinguishable. Second, it was assumed that the 'refined index' of this 
quantum mechanical problem is related to the equivariant volume (i.e. the integral of y^"^^ 
over A^ri) by a simple overall multiplicative renormalization required for consistency with 
angular momentum quantization. Third, it was important that, in the case of loosely bound 
constituents relevant for wall-crossing, the phase space A1„ was compact and the action of J3 
had only isolated fixed points, corresponding to collinear multi-centered solutions along the z- 
axis. Under these circumstances, the localization theorem of Duistermaat and Heckman [22] 
can be used to express the configurational index g[\ai\^ t), or rather its refined generalization 
g{{ai}] y,t), as a finite sum over the fixed points of J3. This result was found to match the 
prediction of the known wall crossing formulae from supergravity [23, 24] and mathematics 
[19, 20] in all cases where it was tested. As we shall show in this work, provided A^.„ is 
compact and the fixed points are isolated, the multiplicative renormalisation postulated in 
[17] can in fact be derived from the Atiyah-Bott Lefschetz fixed point theorem [25, 26, 35, 34], 
a quantum-mechanical version of the Duistermaat-Heckman formula. 

In this paper, we show that the same approach can be used to find the spectrum of 
multi-centered black hole solutions at a generic point in the moduli space, given the indices 
associated to single-centered black holes. There are however some important differences: 

1. The set of collinear multi-centered solutions must not only satisfy the BPS equilibrium 
conditions of [1], but also lead to a regular metric. This condition was automatically 
satisfied for loosely bound states near a wall of marginal stability, but needs to be 
checked when computing the index at a generic point in moduli space and for generic 
charges. This condition is expected to rule out all but finitely many decompositions of 
the total charge 7 into a sum Y2'i=i [23]- It can also rule out certain connected com- 
ponents of M.n even when the decomposition 7 = Yll=i is allowed, see §3.3.1 for an 
example. A necessary (but not sufficient) condition is that the collinear configuration 
be regular along the 2;-axis. 

2. For some range of charges carried by the centers, the phase space includes 'scaling 
solutions', i.e. regions where the relative distances between a subset (or all) the centers 
can become arbitrarily small [2, 27, 28, 23, 29]. As a result, the space A4.n is non- 
compact, and the sum over regular collinear configurations fails to produce a sensible 
answer. In particular, it does not have a finite limit as y — )■ 1, and cannot be interpreted 
as a sum of characters of 5*0(3) (nor of its double cover SU (2)). However, despite being 
non-compact, turns out to have a finite symplectic volume, suggesting that it may 
admit a compactification. In the case of 'dipole halo' configurations, introduced in 
[16, 29], it is straightforward to construct the compactification explicitly. The resulting 
space M.n still admits an hamiltonian action of 5*0(3), but the fixed points of J3 are 
no longer isolated, in particular there is a codimension 4 submanifold of fixed points 

allow for y ^ 1, at least as the intermediate steps, in order for localization methods to apply. It would be 
interesting to understand the dependence of Tr' {—yf"^^ on the string coupling and other hypermultiplet 
fields. 
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where the total angular momentum vanishes, and which parametrizes scaling solutions. 
We shall assume that always admits a compactification A4n, although we shall 
not require the details of its construction. 

3. Due to the fact that the action of J3 on Ain has non-isolated fixed points, the equiv- 
ariant volume and equivariant index are no longer related by a simple multiplicative 
renormalization. While the 'refined index' could still be in principle computed by lo- 
calization using the Atiyah-Bott theorem, this would require a complete understanding 
of the compactification Ain which we have not achieved so far. Instead, we propose a 
'minimal modification hypothesis' which determines the contribution of these scaling 
regions from that of the regular, well-separated coUinear fixed points. Our prescrip- 
tion amounts to requiring that scaling solutions contribute with the smallest possible 
angular momentum compatible with the final result being a character of SU{2). This 
prescription is motivated by the fact that classically (ignoring angular momentum 
quantization), scaling solutions carry zero total angular momentum. While we do not 
have a proof of this hypothesis, we shall demonstrate that it is consistent with wall- 
crossing and with the split attractor flow conjecture [23]. For a special class of 'dipole 
halo' configurations, where the moduli space of multi-centered solutions is fully under- 
stood, we shall verify that our prescription agrees with the geometric quantization of 
Ain performed in [16, 29], for an arbitrary number of centers. 

We shall now summarize our proposal. We denote by Q^^^{a,y) the index Tr'(— y)^"^^ 
carried by a single centered black hole with charge a (here, Tr' denotes the trace after 
factoring out the bosonic and fermionic zero modes). Since single centered black holes carry 
zero angular momentum[30, 31], we expect that fl^^j{a,y) is independent of y. However we 
shall not make use of this information in our analysis, and proceed with general Q^^f{a,y). 
Let us denote by flj-efil,y) = Tr'(— i/)^''^ the total contribution to the index from single and 
multi-centered black hole solutions carrying total charge 7, and by 

-1 y-y^^ 



y,U _ y 

First consider the case when there are no scaling solutions. Then our proposal for firef (7, y) 
is: _ ^ 

f^ref(7,2/)= 5Z ^7771— ry^rcf («!,•• -,««;!/) ^^fef("i;y)- ■■^rrf("n;i/)- (1-2) 

Here F is the charge lattice, Aut({aj}) is the symmetry factor appropriate for Maxwell- 
Boltzmann statistics,^ and VL^^^{a,y) is the 'rational refined index', related to the refined 
index by 

m\a 



^Aut({ai}) is defined as the order of tlie subgroup of the permutation group of n elements which preserves 
the ordered set (ai, . . . , a„), for a fixed (arbitrary) choice of ordering. Thus if the set {ai} consists of ri 

h 



copies of r2 copies of (^2 etc. then |Aut({ai})| = Ylk ^k- 
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The coefficient g^ei is the refined index of the configurational degrees of freedom of n-centered 
BPS black hole solutions. By localization, it evaluates to 



(1.4) 



where aij = {ai,aj) is the (Dirac-Schwinger-Zwanziger) symplectic inner product between 
ttj and aj, the sum over p represents sum over all coUinear solutions to the BPS equilibrium 
conditions (2.24), (2.10), and s{p) takes value ±1 as determined from eq.(2.29). As mentioned 
above, the set of possible decompositions of 7 that contributes to the sum (1.2) is expected 
to be finite [23]. We can recover the integer invariant f2ref from the rational invariant firef 
using the inverse formula [17] 

aef (7, y) = Yl /^(^) ^"'(^ - y'^^y"" - yn-'^rAi/m, y^) (i.s) 

m|7 

where /i(m) is the Mobius function. Expressing this in terms of fi^j using (1.2), (1.3) we 
can arrive at an expression of the form 

firef(7,2/)= E G({A},{m,};y)n^fef(A,2/"^'), (1-6) 

{ft6r},{m,ez} i 

for some function G. This expresses the total index in terms of the indices of single centered 
solutions. Of course, the sum in (1.6) also includes the contribution from single centered 
black holes given by f2^£(7,?/). 

In the presence of scaling solutions, (1.6) cannot be the full answer for the following 
reason. If we follow the procedure outlined above ignoring the presence of scaling solu- 
tions, and denote the corresponding functions G by Gcou, we shall find that the functions 
G'coii({/3j}, {'^i}; y) are not Laurent polynomials in y, i.e. finite linear combinations of 
with integer m. Hence the corresponding expression (1.6) cannot be interpreted as the gen- 
erating function of the spectrum of a quantum mechanical system with quantized angular 
momentum. Our prescription for taking into account the effect of scaling solutions is to 
modify (1.6) to 

^ref(7,y) = G^coll({A},{m,};i/) n (f^fef(A,y"^0 +^scaling(A,Z/"^0) , (1-7) 

{/3ier},{miGZ} i 

where l^scaiing(a, 2/) is given by 

acaiing(a,2/) = ^({A},{m.};2/) Ha^'eflA,?/"'), (1-8) 

{/3i6r},{migZ} i 

for some function H{{f3i}, {mi}; y). To determine H we substitute (1.8) into (1.7) to express 
the latter equation as 

acf(7,y)= G({A},W;y) narf(A,2/"^0, (1-9) 

{|3i6r},{mjeZ} i 
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for some functions G. We fix H by requiring that G{{(5i}, {mj}; y) be a Laurent polynomial 
in y. The ambiguity of adding to H a Laurent polynomial is resolved by using the minimal 
modification hypothesis, which requires that H must be symmetric under y ^ y^^ and 
vanish as ?/ — oo. An iterative procedure for determining the functions H and hence G has 
been described in §5. 

One advantage of our construction compared to the split attractor flow conjecture of 
[23] is that it allows us to compute the contributions from scaling solutions. Since the 
latter are usually stable across walls of marginal stability (as illustrated for three centers 
in §2.3 and Fig.l below), their contributions cannot be obtained from attractor flow trees. 
This advantage is however mitigated by the fact that we do not know how to determine 
a priori which decompositions 7 = lead to regular multi-centered solutions, except 

by checking (2.10) numerically. It should also be emphasized that in cases where the BPS 
spectrum is described by some quiver quantum mechanics with non-zero superpotential, 
the configurational degrees of freedom carried by scaling solutions do not exhaust, by far 
[29], the exponentially large number of states present on the Higgs branch [23]. These Higgs 
branch states have macroscopic entropy and should be included as part of the single-centered 
configurations counted by f2^f(7). 

The rest of the paper is organized as follows. In §2, we review the structure of the 
phase space of multi-centered BPS solutions in A/" = 2 supergravity, establish the finiteness 
of its symplectic volume, and compute the classical and quantum refined index (also known 
as equivariant volume and equivariant index) by localization. In §3, we investigate several 
examples of three-centered solutions in a simple one-modulus supergravity, paying attention 
to the regularity condition (2.10) and to the contributions of scaling solutions. In §4 we 
describe our proposal for the index of multi-centered black hole configuration when there are 
no scaling solutions, and show the consistency of this formula with wall crossing and with 
the split attractor flow conjecture. In §5 we describe our proposal for modifying the result 
of §4 in the presence of scaling solutions. In §6 we prove the validity of our prescription 
in the solvable case of dipole halo configurations. Further technical details are relegated to 
appendices: in appendix A we provide a detailed analysis of certain sign rules which govern 
the contributions of coUinear fixed points. In appendix B we prove that in the absence of 
scaling solutions the right hand side of (1.4) can be expressed as a Laurent polynomial in 
y. Explicit computations of the equivariant volume and index for dipole halos with 4 and 5 
centers can be found in appendix C. 

2. The phase space of mult i- centered configurations 

We begin by reviewing some relevant properties of supersymmetric multi-centered black hole 
solutions in A/" = 2 supergravity. Such solutions fall into the stationary metric ansatz 

ds^ = -e^^ {dt + A)^ + e-^^dr-^ (2.1) 

where the scale function U, the Kaluza-Klein one-form A and the vector multiplet scalars 
t"',a = 1, . . . ,ny depend on the coordinate r on M.^. 
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2.1 Equilibrium and regularity conditions 

Let r denote the charge lattice. Locally in the moduli space of the theory a charge vector 
a G r may be split into its electric and magnetic components (p^,Q'a)- Given two such 
vectors a and ol we can define a symplectic product 

(a, a') = gAp'^ - (2.2) 

For n centers located at rl, . . . , fvi, carrying electromagnetic charges ai, . . . , a„ in the charge 
lattice r, the values of the vector multiplet scalars and of the scale factor V are obtained by 
solving the "attractor equations" [1]. 



-2 e-^('^ Im {e-'^^Y ((^ (f))] = /3 + ^ ^^Z^M ' = ""^^ ^^"^^ 



where y(t") = — e'^/^(X^(t), FA(t)) is the symplectic section afforded by the special geometry 
of the vector multiplet moduli space, /C = — \iii{F^X'^ — F^JK.^)^ and is the central charge 

Z, = (7,F(too)) (2.4) 

associated to the total charge 7 = + - ■ ■ + «„• ^(^oo) denotes the value of Y at infinity. The 
constant vector /3 on the right-hand side of (2.3) is determined in terms of the asymptotic 
values of the moduli at infinity by 

/3 = -2Im[e^^<^F(t„o)] • (2.5) 

The locations are subject to the equilibrium conditions (also known as integrability equa- 
tions) [1] 

n 



Y.^=c.. (2.6) 



3=1 '"^J 

where r^j = |rj — rj|, = (pti^aj), and the real constants 

Ci = 2Im(e-^'^Z,J (2.7) 

depend on the the asymptotic values of the moduh. Since = arg Z^, these constants satisfy 
Yl^=i^i = 0- Finally, A is given by 

n 



+ , (2.8) 



r — rJ — ' r — Tj 

1=1 ' ' 1=1 ' 



where *3 denotes Hodge dual in 3 fiat dimensions. The conditions (2.6) guarantee the exis- 
tence of a Kaluza-Klein connection A such that the above configuration is a supersymmetric 
solution of the equations of motion. 

It follows from (2.3) that the scale factor U is given by evaluating the Bekenstein- 
Hawking entropy function S{'^) on the harmonic function appearing on the right-hand side 
of (2.3) [3], 
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In order for the solution to be physical one must require that the scale factor be everywhere 
real and positive [23]^, 



+ >0' V^"^^'' (2-10) 



i=l 



where is the location of the z-th center. This is necessary to ensure that there exists a 
regular solution to the attractor equations (2.3) at all points in M.^. 

2.2 Symplectic structure and equivariant volume 

Leaving aside the regularity condition (2.10) for now, let us denote by Ain{{c(ij}', {q}) the 
space of solutions {fi, . . . , r„} to the equilibrium conditions (2.6), modulo overall translations 
of the centers. Ain is a (possibly disconnected) 2n — 2-dimensional submanifold of ]R^"'~^\A, 
where A is the locus in ]R^'"~^ where two or more of the centers fi coincide. ]R^"~^\A is 
equipped with the closed two-form^ 

1 e"''^ drf . A drf . r^- 1 ^ ^ ^ 

w = - 2^ aij 1^ 13 = 2 ^ ^ ^"^'^ ' 

i<j i<j 

where 6ij,(j)ij are the polar angles parametrizing the direction of the vector fy with respect 
to a fixed unit vector u, for example u = (0, 0, 1). For generic values of Cj, the restriction of 
oj to Ain is non-degenerate^, and provides A^n with a symplectic structure [16]. 

By construction, the symplectic form u is invariant under S0{3) rotations. The moment 
map associated to infinitesimal rotations is the angular momentum 



J 



^E«^^S = ^E^^^- (2-12) 



i<3 



where the second equality follows by using (2.6). After some further algebraic manipulations, 
the norm of J can be written as [16] 



\ 



1,3 



'^This condition is somewhat too strong, since it rules out the 'empty hole' configurations representing 
e.g. the BPS states which become massless at the conifold point [1]. A more accurate requirement would be 
that regions with S* < are shielded by a shell where the moduli lie at singular conifold-type points, such 
that the regions inside the shell can be replaced by flat regions with constant values of the scalars. To avoid 
this complication we shall restrict to charges which satisfy D{ai) > 0, where D is the quartic polynomial 
such that S = ^/D/tt in the large volume limit. 

*Our normalisation differs by a factor of two from the one used in [16]. This ensures that cj/27r has integer 
periods, see §2.5. 

notable exception occurs when Ci = 0, where the space of solutions admits an exact dilation symmetry 
Vij — )■ ETij, along which the two-form uj becomes degenerate. This case corresponds to multi-centered 
solutions asymptotic to AdS2 x S'^, and will become relevant in the analysis of scaling solutions below. 
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For n > 3 centers, the orbits of the 5*0(3) action are generically 3-dimensional, except on 
the two-dimensional subspace A^n;coii C Ain corresponding to colhnear configurations. Re- 
moving this locus, the quotient of A^n\A^n;coii by the 5*0(3) action is a (2n — 5)-dimensional 
Poisson manifold Ain- Let JT" = j{Ain) C be the range^ spanned by the total angular 
momentum j on the space of configurations Ain- The symplectic leaves of Ain are the hy- 
persurfaces Ain{j) with fixed total angular momentum j G If one so wishes, one may 
parametrize Ain by 2n — 5 relative distances (suitably chosen among the n{n — l)/2 radii 
Vij, and subject to triangle inequahties), and the leaves Ain{j) by their projectivization. 
Using Euler angle coordinates 9,(j),a on 5*0(3), the symplectic form on Ain may then be 
decomposed into 

uj = j sin Ode Ad(t)- dj A (dcr + cos 6'd0) + u , (2.14) 

where u is the symplectic form on the symplectic leaf Ai{j). Eq. (2.14) follows by requir- 
ing the invariance of u under the vector fields dcr,d^, {da — cos 6*9^)/ sin 6', generalizing the 
argument in [16] to an arbitrary number of centers. 

Our goal in this work will be to determine the 'refined index' Tr'(— y)^"^^ of the super- 
symmetric quantum mechanics of n-centered BPS configurations, where J3 is the angular 
momentum operator along the z-axis. We defer to §2.5 a detailed discussion of this quan- 
tum mechanics, and focus for now on the classical version of the refined index, the phase 
space integral [17] 

*,...,.,({«.}; .) ^ e'-'' . (2.15) 

where throughout this paper, we denote 

u = \ny. (2.16) 

Such integrals over a symplectic manifold of the exponential of the moment map of some 
Hamiltonian action are well studied in the mathematical literature under the name of 'equiv- 
ariant volume' (see e.g. [32] for a survey). The convergence of the integral (2.15) will be 
addressed in §2.3. The sign (— l)Ei<j (2.15) is inserted for reasons which will be- 

come clear in §2.5. Leaving aside this sign, the equivariant volume (2.15) is expected to be 
a good approximation to Tr'y^"'^ in the classical limit, where all symplectic products aij are 
scaled to infinity and y ^ 1 (this last condition ensuring that the function y^"^^ varies slowly 
on the phase space). 

Using the description of Ain as 5*0(3) x Ain we can rewrite the integral (2.15) as 

tacai({«.};y) = „ Z " / ^"^^ / sin^d^d^dae^^^™^^ / Co"^-' . 

(27r)" '{n-Sy. Jj Jso{3) JMnij) 

(2.17) 

^In general, is bounded, and consists of a set of intervals in M+. The regularity condition (2.10) may 
rule out certain intervals in JJ. 

''We shall use the symbol J3 to denote both the quantum angular momentum operator as well as the 
classical angular momentum (2.12). The correct interpretation should be clear from the context. 
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Carrying out the angular integral, we arrive at 

^?classical({«.};i/) = (-1)S.<. / dj ^^^^^^^ ({«,}, J ) (2.18) 

where 



J 



^classical({ai}, j) = ,^ s^_^, ^ / ^ (2.19) 

(27r)" 3(^_3)! y^^^^.^ 
is the symplectic volume of the leaf A^„(j). 

2.3 Scaling solutions and finiteness of the equivariant volume 

In order to assess the convergence of the equivariant volume (2.15), it is useful to start with 
the simplest three-body case, which was discussed in detail in [16]. Recall that the solutions 
to the equilibrium conditions 

ai2 , ai3 "12 , «23 "13 "23 

\ =Ci, \ = C2, = C3 = -Ci-C2 (2.20) 

ri2 rxz ri2 ris 
can be parametrized by [16] 

"12 "23 "31 /ooiN 

^2 = , ^^23 = \ , ^"13 = , (2.21) 

P P + C2 P - Ci 

where p runs over the subset of M satisfying r^j > and the triangular inequalities ri2 < 
'"is + ^"23, ''"23 < '''12 + ^"13, ''"13 < ''"12 + ^"23- In general, the allowed range of p consists of 
at most two intervals, possibly reaching ±00, whose finite endpoints correspond to coUinear 
configurations (see Fig. 1 for a pictorial determination of the allowed range of p as a function 
of the constants q). Trading p for the total angular momentum j using (2.13), one finds 
that the range J of the latter is always a bounded interval [j_, j+]- 

When ai2, ^23 and aai are positive and satisfy the same triangular inequalities, the 
region p — )■ 00 is included in the allowed range, and corresponds to scaling solutions where 
the three centers come arbitrarily close to each other [16]. From the second equality in (2.12) 
it follows that such configurations have vanishing total angular momentum j in the strict 
limit p — )■ 00, hence j_ =0. 

Even though the space A^3 is non-compact when such solutions are allowed, its equivari- 
ant volume S'ciassicaill^i}; v) (in particular, its symplectic volume) is in fact finite. To see this, 
note that the reduced symplectic space A^3(j) consists of a single point when j G [j_, j+], 
or is empty otherwise. Thus, 

In the presence of scaling solutions, j„ =0, but the integral is still convergent. In fact the 
space A^3 may be compactified by adding one point , corresponding to scaling solutions in 
the strict p = 00 limit. Since j = in this limit, this point is fixed under the action of 
50(3). 

We now turn to solutions with arbitrary number of centers. In the special case where 
the centers can be ordered via a permutation a of {1, . . . , n} such that ao-{i)o-(j) > whenever 
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Figure 1: Phase diagram for 3-center configurations, in the plane (x, y) = (—§(02+03), ^(02 — 03)). 
The shaded area corresponds to the values of q ri2 reached by varying the position of the third 
center, keeping the position of the centers 1 and 2 and aij fixed. The range of ri2 (and therefore, 
of the parameter p in (2.21)) at fixed values of Ci can be determined by intersecting the shaded 
area with a radial line extending from the origin through the point Cj (Thus if the radial line 
passing through and c meets the boundary of the shaded region at c^^^ and c^^\ then the range 
of ri2 is from |c^^)|/|cl to |c^^^|/|c|. On the other hand if the origin falls inside the shaded region 
then the range of ri2 varies from to |c^^^|/|c|). The blue, red and yellow lines correspond to 
collinear configurations in order 132, 213, 321 (or their reverse), respectively. The lines with 
aig{x + iy) = 7r/2, — 57r/6, — 7r/6 denote the locus where ci = 0,C2 = 0,C3 = 0, respectively, i) 
Q12 > 0, 023 > 0, ai3 > 0: solutions exist in the sector ci > 0, C3 > ii) ai2 > 0, Q23 > 0, ai3 < 0, 
012 > «3i + 023 (non-scaling regime): solutions exist in the sectors ci > or C2 < iii) 012 > 
023 > 0, ai3 < 0, ai2 < a3i + 023 (scaling regime): solutions exist for all values of Cj. The figures 
were produced using (012, 023, 013) = (1, 3, 2), (6, 3, —2), (1, 3, —3) in cases i), ii), iii), respectively. 

i > j, ® it is straightforward to see that the bounded from below by a non-zero rmm! 

therefore scaling solutions do not exist. Away from walls of marginal stability or threshold 
stability^, the r^j's are also bounded from above, and the phase space M.^ is therefore 
compact. In such cases, the equivariant volume (2.15) is manifestly finite. 

If on the other hand, the centers cannot be ordered such that a^-j-j) ^.(j) > whenever 
i > j, then there may exist a subset A of (1, . . . ,n) for which there exist vectors fi G 
{i G A) satisfying 

^^=0, V leA. (2.23) 

*This is equivalent to the condition that the associated quiver (with n nodes and arrows from i to 
j if aij > 0, or from j to i if aij < 0) contains no oriented closed loop. This condition is sufficient for the 
absence of scaling solutions, but not necessary, as illustrated by the 3-center case. It is in particular obeyed 
when all charges ai belong to a positive cone in a two-dimensional charge lattice, the case studied in [17]. 

^Walls of marginal or threshold stability correspond to loci where J^ieA Ci = for a proper subset A of 
{1, . . . , n}, such that the centers in the subset A can move off to infinity. Threshold stability corresponds to 
the special case where the charge 7^ = X^jgA '^i '-'^ subset A and the charge 7 — 7^ of its complement 
are mutually local [33]. 
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In this case the centers in the subset A may form a scahng solution, i.e. may come arbitrarily 
close to each other. In this case, the space Ain is non-compact. However, we shall now outline 
a proof that the equivariant volume (2.15) is in fact finite even in this case. 

We start by considering potential divergences from 'maximally scaling solutions', where 
all centers come arbitrarily close to each other. Since from eq.(2.11) the symplectic form u 
is invariant under rij — )■ evij, the volume element o;""^ is superficially expected to behave 
as de/e in the region e — 0, suggesting a logarithmic divergence of the equivariant volume. 
However, notice that in the e — )■ limit, the symplectic form u becomes degenerate in 
the direction corresponding to overall dilations of the configuration (as noted in footnote 5 
above). This dilation is an exact symmetry of the equilibrium equations (2.6) at Cj = 0, but 
is broken for finite e at non-zero q. Thus, the naive estimate cj"~^ ~ de/e: must be modified 
by an integer power^° k > 1 of the dimensionless parameter ce (where c denotes a function 
of the Cj homogeneous of degree one). In this case the volume element cj"~^ ~ {ce)''de/e 
vanishes at Cj = 0. Thus, the contribution of the region e — )■ to the equivariant volume is 
finite, and there is no divergence coming from 'maximally scaling solutions'. 

We should also consider situations where m of the n centers come arbitrarily close to 
each other. The phase space near this configuration is locally a product of the phase space 
of an [n — m + 1) centered configuration, in which all the m close-by centers are regarded 
as a single center, and that of a m-centered scaling solution. The symplectic volume of the 
former is manifestly finite, while the symplectic volume of the latter is finite by the argument 
given earlier. This can be easily generalized to the cases where there are several subsets, each 
containing a set of centers which come close to each other. Thus, this reasoning indicates 
that the equivariant volume (2.15) is finite even in the presence of scaling solutions. In §6, 
we shall verify this claim in the case of dipole halo configurations, where the phase space 
admits a natural compactification into a toric symplectic manifold, whose volume can be 
evaluated explicitly. It is an important problem to construct a compactification of Ain in 
the general case involving scaling solutions. 

2.4 A fixed point formula for the equivariant volume 

Having established that the equivariant volume (2.15) is finite, we shall now evaluate it by 
localization. For this purpose, we need to determine the fixed points of the action of J3 on 
Ain- In general, the fixed points correspond either to coUinear solutions along the z axis, 
or to the coinciding limit of scaling solutions approaching a point on this axis. First we 
restrict to the case where the charges do not allow for scaling solutions. In this case, Ain 
is compact and all fixed points are isolated. 

CoUinear solutions lying along the z axis satisfy a one-dimensional version of the equi- 

^°In fact, using the decomposition (2.18) of the equivariant volume and the second relation (2.12), one 
finds that the volume element jdj along the angular momentum direction scales as c^ede. If the symplectic 
volume of the reduced phase space Mnij) has a finite non-zero value at j = 0, which is the case for general 
3-ccntcrcd configurations or for dipole halo configurations with an arbitrary number of centers, it follows 
that the vohmrc clement goes as ede near e — 0. 
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librium conditions (2.6), 



E^^=^- E-^ = 0' (2-24) 



I Z-i 



where the second equation fixes the translational zero-mode. These equations then select 
the critical points of the 'superpotential' [17] 

W{\ {^i}) = - ^ sign[2j - In \zj - Zi\- - A/n)^, (2.25) 

as a function of n + 1 variables A, 2i, ■ ■ ■ Zn. In the vicinity of one such fixed point the 
angular momentum J3 and the symplectic form bj take the form 

J3 = - ^ ttij sign[2;j - z^ - -Mij{p) {xiXj + yiVj) H , ^ = 2 ^^'^^^"^ ^^^Vj^ > 

(2.26) 

where Mij{p) is the Hessian matrix of W{\, {zi}) with respect to 2:1, . . . , 2;^, and (xj, Ui) are 
coordinates in the plane transverse to the 2;-axis at the center 2, subject to the condition 
'Ylii = Y^Ui = 0- Except for the overall translational zero-mode, the matrix Mij is non- 
degenerate, and the critical points are isolated. Since Ain is compact, we have a finite set 
of collinear fixed points p. It will be useful to note that det M = — det M, where det M 
denotes the determinant of the (n — 1) x [n — 1) matrix M with the first column and first 
row removed (i.e. the Hessian of 1^ as a function of Z2, . . . , Zn), and M is the Hessian of W 
as a function of the n + 1 variables X, zi, . . . , Zn [17]. 

Since A4n is compact and the fixed points are isolated, the Duistermaat-Heckman for- 
mula [22] allows to express the equivariant volume (2.15) as a sum over fixed points. The 
contribution of each fixed point p is given by the formal Gaussian integral^^ obtained by 
replacing the moment map by its quadratic approximation around p, and the symplectic 
form by its value at the same point. Thus the net contribution is given by 



(_l)E.<,^.+n^i r " (detM(p)) 

^classical ( I «i I ;Z/) = (^^^^^=1 J 1 1 d?/i 



p i=2 



n n 

exp |z/ aij sign[zj - Zi]-'^^ Mij{p) {xiXj + ViVj)^ 



(2.27) 



i<j i<j 



After carrying out the Gaussian integral over Xi,yi, we arrive at [17] 



^7ciassicai({a.}; v) = ^ (^i^^^^^i E '(p) ^^^"^ ''"^'''^'^ ' (2-28) 



where the coefficient s{p) is given by 

s{p) = sign det M(p) = -sign det M{p) . (2.29) 

^^This Gaussian integral is ill-defined since the quadratic form is in general not positive definite, but the 
Duistermaat-Heckman formula guarantees that formal Gaussian integration leads to the correct result. 
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—s{p) is recognized as Morse index of the critical point p of W. It is worthwhile noting that 
the prefactor (2z/)"~^ in the denominator of (2.28) originates from the determinant of the 
quadratic form in the 2n — 2 transverse directions to the fixed point. 

In the presence of scaling solutions, is non-compact and the Duistermaat-Heckman 
formula does not apply directly. In the examples that we shall study in §6, it appears 
that A4n admits a natural compactification Ain, which introduces additional, non-isolated 
fixed points. Their contribution can in principle be evaluated using a suitable generalization 
of the Duistermaat-Heckman formula [34]. For lack of a complete understanding of the 
compactification A^^, we shall not pursue this approach. Instead, one may give a 'minimal 
modification prescription' which determines the contribution of the non-isolated fixed points 
by requiring that the Fourier transform of the equivariant volume J^gciassica.i{{c(i}, i^)e^"^'"dz/ 
becomes a piecewise polynomial function of m with compact support, as must be the case 
for any compact symplectic manifold with a Hamiltonian action [22]. Since our interest lies 
eventually in the equivariant index rather than the equivariant volume, we shall not explain 
this prescription here, referring to §6 and §C for direct computations of the equivariant 
volume for dipole halo configurations, and identifications of the corresponding fixed points. 

Finally, we return to the regularity condition (2.10). In the case considered in [17], 
namely multi-centered solutions whose charges all lie in a positive cone of a two-dimensional 
lattice, it appears that this condition is automatically satisfied near the wall. In general, this 
need not be so. On physical grounds, we do not expect that this condition should cut out 
regions of a given connected component of Ain- On the other hand, it is quite possible that 
this condition could forbid certain connected components of Ain altogether. In this case, 
we simply need to omit the corresponding fixed points from the sum (2.28). We shall see an 
example of this phenomenon in §3.2 below. 

2.5 A fixed point formula for the refined index 

We now turn to the problem of quantizing the configurational degrees of freedom of n- 
centered black hole solutions and the computation of the refined index Tr '(— y)^"^^, of which 
(2.28) is supposed to be the classical approximation. 

We start by recalling and expanding upon some general aspects of the quantization of 
BPS multi-centered black hole solutions, first laid out in [2, 16]. As in these references, we 
assume that the Hilbert space of BPS states, in principle obtained by first quantizing the 
full configuration space of A/" = 2 supergravity (or more generally, string theory) and then 
restricting to the subspace annihilated by 4 supercharges, can be alternatively obtained by 
first restricting to classical solutions annihilated by the same number of supercharges, and 
then quantizing. 

Since each BPS solution of the type (2.3) preserves 4 supercharges, one expects that the 
locations should be quantized into operators in a quantum mechanics with Af = 4 super- 
symmetries, such that the classical supersymmetric vacua are in one-to-one correspondence 
with solutions to the equilibrium conditions (2.6). The final result should not be sensitive 
to the details of the potential in this quantum mechanics as long as the space of classical 
supersymmetric vacua is given by (2.6). One way to obtain the Lagrangian is to consider 
the Coulomb branch of an A/" = 4 supersymmetric quiver quantum mechanics with n vec- 
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tor multiplets (each including 3 real scalars r^, a U{1) gauge field Ai and four fermions ipi) 
and, for each pair of centers i,j, \aij\ chiral multiplets with charge sign(ajj) under Ai and 
— sign(ajj) under Aj. On the Coulomb branch^^, the chiral multiplets are massive and can 
be integrated out at one-loop. The result is an A/" = 4 supersymmetric Lagrangian for the 
positions and their fermionic partners ipi, with a potential which vanishes when the 3n 
scalars satisfy the equilibrium conditions (2.6). In addition, there is a first order coupling 
J A to the one-form [2] 

^ ^ «ii (1 - cos 9ij) d(f)ij , (2.30) 

i<3 

where u is the same unit vector used in (2.11). The degrees of freedom associated to the 
center of mass motion are decoupled and can be removed by setting J2i = J2i ''Pi — 0' ^o 
that the bosonic configuration space consists of (M^"'~^\A), where A is the coinciding locus. 
Quantum mechanically, the Hilbert space consists of square integrable sections oi S® C over 
M^"~^\A, where 5* is the trivial bundle of rank 2^"~^ over M^""^ obtaining by quantizing 
the 4(ra — 1) fermionic modes and £ is a complex line bundle with first Chern class^^ 
Ci(£) = uj = dA. BPS states correspond to supersymmetric ground states of this A/" = 4 
quantum mechanics. 

For the purposes of studying these ground states, one may return to the classical 
Coulomb branch Lagrangian, integrate out the fiuctuations transverse to the zero-energy 
submanifold M.n C M^""^, and then quantize the system. While we have not studied this 
problem in detail, we assume, in the spirit of [16], that BPS states correspond to harmonic 
spinors on M.^ i-e. sections of annihilated by where 5* is the total spin bundle^^ of 

A^n, £ is the restriction of C to M.n and D is the Dirac operator for the metric g^y induced 
from the fiat metric on ]R^""^\A, twisted by the line bundle C The Dirac operator D decom- 
poses a.s D = + D_ , where D± maps S± to S'=p where S = S+ (B is the decomposition 
into spinors with positive and negative chirality. Moreover, the action of 5*0(3) on lifts 
to an action of SU{2) on S± CB3 C. The refined index Tr'(— y)^"'^ is then given by 

g,,i{{ai}-y) = TiKcrDA-yf'' + KcroA-y?'' ■ (2.31) 
We shall assume that Ker_D_ = 0, so that the refined index gj.cf{{ai}; y) reduces to 

^7rcf({a.};y) = TrKcrD+(-y)'-'' - TrKerD_(-y)'^' • (2.32) 

^^The quiver quantum mechanics also has a Higgs branch. In some cases, including the wall-crossing 
problem considered in [17], the BPS spectra on the Coulomb branch and Higgs branch are identical, but this 
is not so in general, and the Higgs branch states cannot always be viewed as multi-centered black holes. 

-^■^The normahsation chosen in (2.11) ensures that uj e iJ^(M'^"^'^\A, Z), moreover M'^"^-^\A is simply 
connected so C is uniquely defined. The unit vector u in (2.30) correspond to a choice of Dirac strings for 
the line bundle L. 

^^We assume that A^„ admits a spin structure, unlike a generic symplectic manifold which only admits a 
spiuc structure [35]. By integrating out the massive fermions around the supersymmetric vacua, it should 
be possible to construct the spin bundle S* as a rank 2"~^ subbundle of S, but we have not investigated this 
in detail. 



1^ e'^^'^rf^d^M^ 



u 
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When Ain is Kahler, as is the case for = 2, 3 or for the dipole halo configurations 
studied in §6, and Oij large enough, the fact that KerD_ = can be proved as follows. 
In this case S+ and S_ are isomorphic to A^™^(T(J,°'^^A^„) » K^^^ and A°<^'^(T^°'^^A^„) O 
K^^"^, respectively, where K is the canonical bundle of Ain, i-e. the complex line bundle 
of (n — 1,0) forms (the square root exists since A^n is assumed to be spin). In this case 
the Kodaira vanishing theorem states that the cohomology groups H^'^''^\Ain, ^ ® K^^'^) all 
vanish except for q = 0. This shows that KerD_ = 0. In general, we do not how to prove this 
assumption, but it is supported by the fact that it leads to results in agreement with wall- 
crossing. Since the refined index Tr '(— y)^"'^' is not protected^^ in A/" = 2 supergravity, for 
the purposes of computing the index Tr'(— l)^"'^ it suffices to make the weaker assumption 
that TrKcrD_(-l)'^=^ = 0. 

Under this assumption, 5'rcf({"j}; —y) reduces to the equivariant index of the Dirac op- 
erator D. Using the Lefschetz fixed point formula established in [25, 26, 38, 34]^^ , assuming 
that Ain is compact, the equivariant index can be written as the integral 



g,,ii{ai};-y)= / Ch(£, z/) A(7W„, z/)U_i (2.33) 

where Ch(£, u) is the equivariant Chern character of the line bundle C, and A(A^„, u) is the 
equivariant A-genus of A^n, defined by 

Ch(£, .) = exp + . i(A1„. .) = det ( .J^y^ff^^) ) ■ (2-34) 

Here, J3 is the moment map of the action of rotations around the z axis, L is the endomor- 
phism of the holomorphic tangent bundle T^-^'^'^Ain induced by the same action, TZ is the 
curvature two-form on T'^^''^)Al„,and A\p denotes the degree 2p part of a multi-form A. The 
integral (2.33) can be evaluated by localization [38], leading to 

grediai}; -y) = / — T7,,r/\r a.< fixed ^ — (2.35) 
infixed Eu(iVAIIi^^d) p 

where Aln™*^^ C A^„ denotes the fixed point locus, of dimension 2p, and 

<fixcd\ ■ I ~ ^ 1 



Eu(Ar7Wf ) = det \^2uL + — 7^J (2.36) 

is the equivariant Euler character of the normal bundle of A^^^°*^. 

In the absence of scaling solutions, Atn is compact and all fixed points are isolated, so 
NAi^'^'^ = TAtn and the Euler character cancels the numerator in the A-genus. Moreover, 

^^In N — 2 gauge theories at low energies, it may be possible to compute the protected spin character 
(PSC) [36] by quantizing the space of classical multi-centered Abelian dyons along the lines of [37]. Since 
the PSC is protected, it is plausible that the configurational index will be directly equal to the equivariant 
index (2.32), without the need to assume that KerZ?^ vanishes. The vanishing of KcrZ?_ = in this context 
is closely related to the weak positivity conjecture of [36]. 

^^We are grateful to M. Vergnc for guidance into the math literature. 
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Js = I J2i<j Q^ijSign[zj — Zi] and the operator L, representing the action of J3 on T(^'°^(A^„), 
has eigenvalues ±1, with detL = s{p). This leads to the following explicit formula for the 
quantum refined index: 

(7ref ({«.}; -y) ={y - 1/2/)""' E ^(P) 2/^-^ "'^^'^°'^^-^'^^ ■ (2.37) 

p 

Hence, after changing^® y — )■ —y, 

^ref ({«.}; y) = ^~^^^'^'^^p ' Yl ^(P) ^^'"^ '^^^^'^^'^-'^ . (2.38) 



In particular, the refined index (2.38) is related to the equivariant volume (2.28) by an overall 
rescaling by a factor of {u/ sinhi^)""^. This multiplicative renormalization was postulated 
in [17] on the basis of angular momentum quantization. As the present derivation shows, 
this multiplicative renormalization in fact follows from the Atiyah-Bott Lefschetz fixed point 
formula for the equivariant index of the Dirac operator on Ain- 

Eq. (2.38) will be our main tool for computing the index of multi-centered black holes 
in the absence of scaling solutions. Although it was derived under various assumptions, we 
shall carry out various consistency checks which build our confidence in this formula. In the 
presence of scaling solutions, the space Ain is non-compact and its compactification includes 
non-isolated fixed points for the action of J3. Our limited understanding of the geometry of 
the fixed submanifold Ai^'^'^ prevents us from computing the index using (2.35), nevertheless 
in §5 we shall give a prescription for computing the contributions from these non-isolated 
fixed points, from the knowledge of the isolated ones. 



3. Case studies in a one-modulus supergravity model 

In this section, we analyze several examples of three-centered solutions in the context of a 
simple M = 2 supergravity model with a single modulus. This model, introduced in §3.1, 
arises by compactifying type IIA string theory on a Calabi-Yau 3-fold X with 62 = 1 in the 
large volume limit. The first example in §3.2 illustrates the importance of the regularity 
condition (2.10), while the second example in §3.3.1 shows the necessity of including contri- 
butions from fixed points associated with scaling solutions. This example is a particular case 
of a general class of dipole halo configurations which will be discussed further in §6. The 
example in §3.3.2 illustrates the role of the regularity condition (2.10) in deciding whether 
a given scaling solution is physical or not. 

3.1 One-modulus model 

We consider a simple supergravity model with one vector multiplet, governed by the prepo- 
tential 

F(X°, X^) = -K {Xy/6X^ + ^A{X^y + ^i3X°Xi . (3.1) 

-'^''The same localization techniques show that the index of the untwisted Dirac operator Jj^ A = 
{y — l/y)"~^ J2p ^{p) vanishes, in particular is integer, consistently with the fact that A4n is spin. 

^^The fact that y — —y changes the refined index by an overall sign shows that for fixed aij and n, all 
states carry the same parity of 2J3. 
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This model arises in the large volume limit of type IIA string theory compactified on a 
Calabi-Yau 3- fold X with 62 = 1- In this case n = J^oo^ is the cubic self-intersection of the 
unique generator u of of H^{X, Z), B = / a; A C2 is an integer such that B = —2k mod 12 
and A is an half- integer constant such that A = k/2 mod 1.^^ The ratio X^/X^ is the 
complexified Kahler modulus t: 

t = B + iJ=^, (3.2) 

with B the NS 2-form potential and J the real Kahler modulus. The central charge of a 
BPS state with electromagnetic charges {p^,qA) is given by 



= e^/^X^ 



^Kt^p^ - l-Kt'^p^ - tqi - go 
2 



(3.3) 



where e'^ = i^X^Fa - X^Fa) = |X°|2 is the Kahler potential, and 



B ^ B 

1i = (li- - go = go - • (3.4) 

With the above conditions on A^B, the charges p^, qi and go are quantized in integer 
units[39] and represent (up to a sign) DQ, DA, D2 and DO-brane charges. We shall denote 

a = {p^,p^,qi,qo) , a = (p°, p\ gi, go) . (3.5) 

Note that given a pair of vectors a, a' we have {a, a') = {a, a'). Solving the attractor 
equations gives for the Bekenstein-Hawking entropy tt | Za (Xfttr ) P of a large black hole in 
this model [40]: 



5 = 7rv//^(pO,pi,g~i,go), (3.6) 
where D is the quartic polynomial 

•^,^1^2 ^0^. ^1 ;^2(^0\2 r^US^^ ^0/;;\3 



^2 

L)(p°,p\gi,go) = — 



kP' kP' Kj'^ Kj K/^ 



(3.7) 



For a multi-centered black hole solution with charges the consistency condition (2.10) 
can now be expressed as 



^ + >0, V f6M^ (3.8) 




For single centered supersymmetric black holes D(p°,p^, gi, go) > implies (3.8) for J ^ 0. 
For multi-center black holes (3.8) does not follow from regularity of the near-horizon regions 
but must be checked independently. 

We shall conclude this subsection with some comments on the effect of a' corrections 
to the supergravity action. One source of corrections originates from world-sheet instantons 



These congruence conditions are special cases of the conditions \nabcp''p'^ — Aatp'' ~ mod 1, 
\i^abcP""p'^'P'^ + ffP" raod 1 for all integer vector see [39]. For the quintic, we have k = 5, A = —11/2, 
B = 50. 



- 18 - 



which contribute to the prepotential F. It is straightforward in principle to incorporate such 
corrections in the analysis of multi-centered solutions. In addition, there are also higher 
derivative corrections to the four- dimensional effective action which greatly complicate the 
construction of exact multi-centered solutions [11, 41]. While these corrections are essential 
in computing the index Q^{a) of certain single centered black holes [12], their effect on the 
functions (^ref I/) governing the index of multi-centered black holes is mild, since the 
detailed information about the action and the solution is needed only to determine whether 
a given coUinear solution has a non-singular metric or not; but given such a solution the 
contribution to 5'rcf({tti}; y) given in (2.38), and its generalization to scaling solutions which 
we shall discuss later, is independent of the action. Thus in our analysis we shall ignore the 
effect of these higher derivative corrections. For explicit computation we shall also set A 
and B to zero for convenience. 

3.2 Non-scaling solutions: D6 — D6 — D6 

We now consider a 3-centered configuration of D6-branes with fluxes, carrying charges e^'^, 
e^'^ and -e(^+^)'^ where U, V are positive integers. This example was studied in detail in 
[23], section 5.2.2. To compute the corresponding charge vectors, recall that the central 
charge associated with a D6-brane carrying charge e^'^ is proportional to 

(3.9) 

Comparing this with (3.3), and using the fact that J C2 A u = B, JuAujAu = k,, we see 
that the first center carries charges 

a,= (^l,U,-Lu'-j^B^-KU' + ^BU^ , (3.10) 

where we have determined the overall normalization of ai by using the fact that for a single 
D6-brane = 1. Similarly the second and the third centers carry charges 





e^- (l H 








24/J 



(3.11) 



a.= {l.V.-\nV^-l-^B^-.V^^l-^BV 

«3 ={-l,-{U + V), \k {U + Vf + Y^B,-^-n {U + - ^^B{U + V) 

so that the total charge is given by 

7 = ai + 52 + as = (^1, 0,f^UV- ^B, -\ UV{U + V)^ . (3.12) 

It is easy to see that for integer U and V, the charges qi are not integers in general, but qt 
computed via (3.4) are. Using aij = aij, the integer symplectic products are then given by 

a2i = lf^{V-U)' + ^B{V-U) , a23 = ltiU' + ^BU , = + ^BV . (3.13) 

12 12 12 

Our goal is to find the index of supersymmetric bound states in the large volume limit. 
We choose concrete example 

K = 6, B = 0, U=l, V = 2, B + iJ = 3i. (3.14) 
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This choice differ from that of [23] in the choice of n, but this is a simple normahzation 
factor and still allows us to compare our findings with the results of [23]. The value of J is 
chosen to be large enough to lie in the large volume chamber. The precise value of J affects 
the numerical values of the solutions to eq.(2.24), but is otherwise irrelevant. We use the 
projective symmetry of special geometry to fix X° = 1. For the values of charges and moduli 
corresponding to (3.14), this leads to 



"13 = 8 



"12 = -1 , "23 = 1 

ci = JIVSTS, C2 = 1^^318 , /3 ^ „ _ _ 
Numerical analysis of (2.24) leads to the following collinear solutions:^'' 



53 ( J_ 2_ 63 _18^ 
6 U59' 53' 53' 53' 



(3.15) 



cr 
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1 




321 
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1.02 


-5.14 


-1 
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321 


-1.02 


-0.195 





-5.14 


-1 


X 


213 


-0.182 





0.974 


-5.49 


1 


X 


312 


-0.974 





0.182 


-5.49 


1 


X 



(3.16) 



In this table, we have displayed the permutation a of 123 specifying the order of the Zi^s for 
a given solution according to Zo-(i) < z„(^2) < z<t(3), the location of the centers 2;o-(i), the value 
of the parameter p in the parametrization given in (2.21), the value of the sign s{p), and in 
the last column, we indicated by a/ solutions obeying the regularity condition (2.10), and 
by X those which do not (it turns out that the region where the discriminant D becomes 
negative intersects the 2;-axis, and therefore can be found by plotting D as a function of z). 
Note that the parameter p takes the same value for two configurations related by a reversal 
of the z-axis. More generally, solutions of (2.6) satisfying the triangular inequalities arise 
from the disconnected intervals —.422 < p < —.363 and —5.49 < p < —5.14 along the p axis, 
and correspond to general non-coUinear solutions with angular momentum in the intervals 
0.444 < < 0.694 and 0.25 < < 0.444, respectively. Only the first interval leads to 
solutions which satisfy the regularity solution (2.10). 

According to our prescription only the solutions marked by ^/ must be included in the 
sum in (2.38), and we get 



5'rcf(ai,a2,a3;i/) 



-1) 



021+«23+ai3 



sinh u 



sinh(a23i^) sinh((ai3-Q;2i)z/) , u = \ny, (3.17) 



in agreement with [23]. Had we ignored the constraint (3.8) and included contributions from 
all the solutions given in (3.16), we would have got 



,l)a2i+a23+«i3 — 3_ sinh(ai3Z/) sinh((a23 + 021)1^) 
sinh u 



(3.18) 



20 



We have presented the solution in the zi = gauge instead of Zj — gauge. 
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This last expression is indeed the result produced by quiver quantum mechanics, but as 
emphasized in [23], quiver quantum mechanics in this case fails to give the correct result. We 
see that the difference between the correct result and quiver quantum mechanics is accounted 
for by the additional constraint (3.8) that must be imposed on the solutions besides (2.24). 

3.3 Scaling solutions 

Consider a three centered black hole solution with the centers carrying charges ai, a2 and 
as. Suppose further that a\2 = 3, q;23 = 4 and a^i = 5. In this case ai2, 0^23 and a^i satisfy 
a triangle inequality, showing that the scaling solution can exist. Thus this provides us with 
a laboratory for studying the role of scaling solutions in the computation of the index of 
multi-centered black hole solutions. 

For now we focus on the constraints imposed by eqs.(2.24) without worrying about 
(2.10); we shall return to this later. Since we have two moduli J and B at our disposal 
we can adjust them to set ci and C2 as we like (within an appropriate range), C3 is then 
fixed to be — ci — C2. Let us use this freedom to choose a point in the moduli space where 
Ci = —-^^jCtij for some positive constant A. We can now look for solutions to (2.24). 
Numerical analysis shows that there are only two coUinear configurations which contribute, 
corresponding to the alignments 123 and 321. Furthermore both contribute with positive 
sign. As a result the net contribution is given by {y — y^^)^"^ (^yCH2+ai3+a23 _|_ y-ai2-ai3-023^_ 

This causes a puzzle since this does not have a finite y — t- 1 limit. We must however 
recall that there are also scaling solutions to Denef's constraints which correspond to all 
three points approaching each other. Clearly this is also a fixed point of J3 if we choose the 
point of approach to lie on the z axis. These however will not show up in the numerical 
determination of the fixed points which assumes from the beginning that the centers have 
finite separation. 

The fixed point associated with the scaling solution has J3 = and hence it contributes 
a constant to the index. There may be additional factors from the integration measure near 
the scaling solution. These can in principle be determined from a detailed analysis of the 
scaling solution. However we can try to guess the contribution from the scaling solution by 
requiring that the total contribution has a finite y — ?■ 1 limit. Our first guess would be 

(7rcf(«l,«2,«3;l/) = (Z/-y~')~' (y--+"13+«23+^-ai2-«13-«23_2) . (3.19) 

For + «i3 + ^23 even, the numerator has a factor of {y — y^^Y, and after canceling this 
against the denominator we are left with a Laurent polynomial in y. This is a sensible result, 
with the coefficient of y"^ counting (—1)™ times the number of states with J3 = m/2. But 
(3.19) does not lead to a sensible spectrum for odd values of ai2 + 013 + 023, since the factor 
of {y — y^^Y denominator is not cancelled. Our proposal is to replace the subtraction 

constant 2 hy y + y^^ in this case so that (y — y~^)~^ factor in the denominator is cancelled 
and we again get a Laurent polynomial in y. 

The above analysis shows that we can predict the existence of and the contribution from 
a scaling solution from the results on the non-scaling solutions which are simpler to find. 
This is the general procedure we shall follow, i.e. assume that the role of the fixed points 
associated with the scaling solution is to essentially make the final result satisfy desirable 
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properties. This by itself will not fix the contribution uniquely. For example in (3.19) we 
could have taken the subtraction term to be + instead of 2. To fix this ambiguity 
we shall make a further assumption, that the extra contribution due to the scaling solutions 
vanish as y — )■ oo. This fixes the correction terms uniquely. We shall call this the 'minimal 
modification hypothesis'. We shall describe the general rule for such replacements in §5. 
While we have no a priori justification for this assumption, it seems to work in all known 
cases. 

For specific choices of the prepotential we also need to verify that the solutions we 
consider satisfy the constraint (3.8). This can be easily implemented by testing (3.8) for each 
of the two coUinear solutions described above and if we find that (3.8) fails for these solutions 
then we would conclude that the subtraction terms associated with the scaling solutions must 
also vanish. Thus there will be no contribution to the index from this configuration. For 
implementing this condition we need to work with specific charge vectors. We shall now 
discuss two examples, - one where this condition is satisfied, and another where it fails. 

3.3.1 Example 1: M - - DO 

Let us take a DQ-DQ-hiane pair with fiux as in the last subsection and a pure DO-charge. 
In particular we choose 

= {l^-U,-KUy2-B/2A,-KUy6-BU/24:), 
a2 = (-1, -f/, K + -B/24, -k U^Q - BU/24), (3.20) 
a-s = (0,0, 0,n), 

for some integer U. In this case we have 

ai2 = 4:KU''^/3 + BU/6, 023 = a^i = n . (3.21) 

We define a = 023 = asi, and assume that ai2 > 0,a > 0. Scaling solutions are expected 
when the triangular inequalities are satisfied, i.e. when 2a > a^- 
We choose the following explicit values of the parameters: 

K = 6, B = 0, f/=l, t = B + iJ = l + 3i. (3.22) 

As example of a system which does not allow for scaling solutions we take n = 3, while we 
take n = 6 for an example of a system that does. For those values, we have 

n = 3: ci = ^y^^, C2 = -^^J~^, (3.23) 
n = 6: c, = f^, C2 = -i^. (3.24) 

For n = 3 one finds collinear fixed points for p = 3.41 with alignment 312 (or its reverse) and 
s{p) = 1, and for p = 11.37 with alignment 132 (or its reverse) and s{p) = —1. A detailed 
study of the full supergravity solutions shows that they satisfy the regularity condition (3.8) 
everywhere outside the centers. The contribution of these fixed points gives 

+ y^^ ~ ^ 

9vci{{ai]\y) = . _ ' (3-25) 
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which is finite for y — )■ 1 as in §3.2. 

For n = 6 one finds a single value of p = 3.99 corresponding to collinear solutions with 
alignment 312 and s{p) = 1 (and its reverse). We have checked that this solutions satisfies 
(3.8) everywhere. Adding the contributions of these fixed points gives {y — y~^)~'^{y^ + y~^), 
which does not have a finite y — )■ 1 limit. Following the procedure described earlier we now 
take the total contribution to be 

(7,ef ({«.}; y) = iy- r ')"' {y' + y~'-2) , (3.26) 

which does give a sensible result for the spectrum. 

We close this subsection with a comment on the structure of the solution in the regime 
< B 1, J ^ 1. For this discussion we parametrize ci and C2 by ci = ^ — rj and 
C2 = — /i — Tj. In the limit described above, > and < ^ 1. For 2a < ai2, there 
are two collinear configurations (up to reversal of the z axis) with alignments 132 and 213, 
parametrized by 



Pirn = + 0{v') and p(2i3)=/i+A/— + 0(r/) , (3.27) 

' ' ai2 - 2a V "12 

respectively. For 2a > ai2, only the second solution leads to physical solutions. In the 
?7 — )■ limit the first solution in (3.27) has a smooth limit, with the center 3 sitting at 
the middle point between the centers 1 and 2, separated by ri2 ~ (ai2 — 2a) /fi. For the 
second solution, we have instead r23 ^ ri3 = y'^^^ + 0(?7°) with ri2 — ai2/p- Therefore 
the point 3 moves off to infinity along the z axis, similarly to what happens at a wall of 
marginal stability. However, since {ai + a2, as) = 0, the index does not change across the 
locus rj = 0. This is an example of a wall of threshold stability [33]. While a non-compact 
direction opens up in A^a at 77 = 0, the equivariant volume (2.28) stays unchanged. In §6 we 
shall see that the locus 77 = (or its analogue for more centers) is very convenient to analyze 
the phase space, even though some of the collinear fixed points sit at infinity. Such 'infrared 
divergent' configurations should however not be confused with the 'ultraviolet divergent' 
scaling solutions, which cannot be removed by moving away from 77 = 0. 

3.3.2 Example 2: Unphysical scaling solutions 

In the example in §3.3.1 the entropy associated with individual centers vanish as can be 
easily seen using eq.(3.6) and (3.20). We now consider another example where each center 
describes a regular black hole with finite entropy. Again for simplicity we take k, = 6 and 
B = and choose 

ai = (0,7,48,0), a2 = (0,9,42,0), ag = (0, -8, -48, 6) . (3.28) 
In this case we have 

ai2 = 138, a23 = 96, a-si = 48 . (3.29) 

Thus the triangle inequality is satisfied and we can look for scaling solution. Taking t = 
B + iJ = 3i, we find that collinear solutions to (2.24) exist for the alignment 132 (or its 
reverse), with s{p) = 1 and p = 62.68. Naively, to the contributions from these collinear 
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fixed points we should add the contribution —2 from scahng solutions, so as to obtain an 
admissible SU{2) character. Numerical analysis shows however that the coUinear solutions 
fail to satisfy (3.8) in some region of the z-axis and hence are not valid solutions. Thus 
the scaling solutions must also be absent. This can be verified explicitly by examining the 
full solution space parametrized by (2.21): in this case the values of p consistent with the 
triangular inequalities range from p = 62.68 (corresponding the above collinear solution) to 
p = oc (corresponding to the scaling regime). We find that throughout this range, the left 
hand side of (3.8) fails to be positive somewhere in the three dimensional space. Thus, the 
regularity condition (3.8) rules out the entire phase space A^s in this case. 



4. Index from non-scaling multi-centered configurations 

In A/" = 2 supersymmetric string theory, typically the index in the sector of charge 7 re- 
ceives contributions both from single centered black holes carrying charge 7, and from multi- 
centered black hole solutions, with individual centers carrying charges ai, . . . , such that 
Y2 - ai = 7. In this section, we shall consider contributions from those multi-centered config- 
urations which do not allow for scaling solutions, i.e. solutions where three or more centers 
can come arbitrarily close. As explained in §2, this requires that there is no subset A of 
{1, . . . ,n} for which we can find vectors fi {i G A) satisfying (2.23). We work at some fixed 
point in the moduli space, and denote by Za the central charge associated with the charge 
a. Using the same logic as in [17] we arrive at the following expression for the (refined) 
index f2i.ef(7; y) = Tr '(— y)^"'^' from multi-centered black hole solutions carrying total charge 
7 (assuming that 7 is a primitive vector of the charge lattice): 

A ./r ,^ 9rei{ai,---,an]y)^refi(^l'^y)---^vei{(^n]y)- (4.1) 

^ Aut({ai}) 

Here, Aut({aj}) is the symmetry factor appropriate for Maxwell-Boltzmann statistics (see 
footnote 2), while 5'ref(tti, • • • ,C(n] y) given in (2.38) is the contribution to Tr'(— y)^''^ from 
n-centered black hole configuration with the centers regarded as distinguishable particles 
with unit index, f^f^f (a; y) is the 'rational refined index', related to the refined index by 

m\a 

It is worthwhile noting that (4.1) and (4.2) exhibit a manifest 'charge conservation property', 
whereby each power of f^^j (a, y) carries charge a and each power of f2ref (a, y"^) carries charge 
ma. The usual protected index is obtained by taking the y — ?■ 1 limit of (4.1), at the cost of 
obscuring manifest 'charge conservation'. If 7 is not primitive, then (4.1) represents 0i.ef(7, y) 
defined in (1.1) rather than r2ref(7,2/)- 

Using (4.1), (4.2) and the relationship between l^ref and VL^ai we can express r2ref(7,2/) 

as 

acf(7,y)= G({A},W;y) n^fcf(A,2/"^0, (4.3) 

{|3i6r},{mjeZ} i 
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for some function G. The i-th term in the sum represents a contribution from configurations 
with a total of centers, with rrii centers of charge Pi. If two or more of the /3i's 

are identical, the total number of centers carrying a given charge Pi is the sum of the 
corresponding m^'s. Finally the sum also contains the contribution from a single centered 
black hole of charge 7, represented by the term Q^^f{'y,y). In order for the right hand side 
of (4.3) to be a bona fide SU{2) character whenever the f2^f(/3i, |/)'s are, the functions G 
which appear in (4.3) must be Laurent polynomials in y. The reader is referred to appendix 
B for a proof of this property in the absence of scaling solutions. 

Eq.(4.1) (or equivalently (4.3)) gives the net contribution to the total index from all 
possible single and multi-centered solutions carrying a fixed charge 7. However it is also 
useful to identify which terms in (4.1) represent the contribution from a specific multi- 
centered solution. We shall assume for simplicity that 7 is primitive so that VL^ciiliU) = 
^rciiliV)- When all a^'s are different then the summand in (4.1) with fif^j replaced by Q^^^ 
represents the contribution to the index from multi-centered black holes carrying charges 
, a„. However when some of the aj's are equal (say ri copies of Pi, r2 copies of P2 
etc.), then there is additional contribution to this index due to the fact that the fi^f (a/m, y™) 
term in (4.2) represents the contribution from m identical centers, each of charge a/m. Thus, 
the contribution to the index from a multi-centered configuration with ri centers of charge 
Pi, r2 centers of charge P2 etc. is given by 

r '^fe r 1 -1 ^ " 

w>im4'"'>i} ^ ^-=' ~y 

nj. (a)_ 



where 



grefiai, ...,an]y) = ^^^(-j^.j^ fi'ref(ai, ...,an]y). (4.5) 



After summing (4.4) over all possible choices of {/3fc}, {r^} satisfying Ylik ^kPk = 7, we recover 
(4.1). As explained in §2, the index (7rcf(«i, • • • , Cin', y) can be computed by localization. Since 
we assume that the phase space contains no scaling solutions, all fixed points are isolated 
and the quantum index is given by (2.38). 

In the remainder of this section we shall carry out various consistency tests of our 
proposal for the index associated with the multi-centered black hole solutions. 



4.1 Consistency with wall crossing 

Let us examine whether our proposal is compatible with wall crossing. For this purpose, 
assume that the moduli are chosen near a wall of marginal stability where the state carrying 

total charge 7 = ai H |-«n becomes marginally unstable against decay into states carrying 

charges 7^1 = XlieA = J2ieB where A and B are two complementary subsets of 

{1, . . . ,n}. At the wall, the phases of the central charges Z^^ and align, and we have 

^c, = -5^c,^0. (4.6) 
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We shall assume that the moduli are chosen on the side of the wall where 

(7A,7s)$^Q >0. (4.7) 

In this region of the moduli space, a class of solutions to (2.24) can be constructed by joining 
two solutions involving the centers in the set A and those in the set B as follows. Let 
and be the cardinality of the sets A and B. Now, choose the relative distances Zi — Zj 
for i,j & A according to a particular coUinear solution involving the charges ai^A, and the 
relative distances Zi — Zj for i,j^B according to the particular coUinear solution involving 
the charges a^gB- Finally, choose the relative separation between the Zi's in the set A and 
the Zi^s in the set B such that 

I I = Ci , 'yAB = (7A, 7b) , (4.8) 

where za = (Sjgyi ^O/'^a zb = {J2i&B ^i)/''^B are the average positions of the centers 
in the sets A and B respectively. Near the wall of marginal stability YlieA 0' ^"^^ 
separation \za — zb\ becomes large. We claim that this configuration satisfies eq. (2.24) in 
the limit (4.6), and can be systematically corrected to an exact solution of eq. (2.24) in the 
vicinity of the wall. 

To see this, we note that if i belong to the set A then eq.(2.24) receives significant 
contribution only when j also belongs to the set A since the distance \za — zb\ computed 
from (4.8) is large in the limit (4.6). Thus the equations reduce approximately to the 
equations for coUinear multi-centered solutions involving the set A only. By assumption 
our solution satisfies the latter equations. A similar argument holds when i belongs to the 
set B. There is however a small caveat stemming from the fact that the solutions in the 
set A are labelled by ra^ — 1 relative distances while solutions in the set B are labelled by 
iriB — l) = (n — nA — l) relative distances. This gives {n — 2) parameters, but (2.24) contains 
(n — 1) independent equations (the sum of the equations over aU i being trivially satisfied). 
Since we cannot adjust (n — 2) parameters to solve (n — 1) independent equations we must 
have missed some equation that determines the relative distance between the points in the 
set A and those in the set B. To find the missing equation we sum (2.24) over alH G A to 
obtain 

ieA J=i ' ^ *' ieA 

Dividing the sum over j on the left hand side into those for which j E A and those for which 
j G B, and noting that the first term vanishes by z o j symmetry, we get 

ieA jeB ' ^ *' i£A 

Approximating zj — Zi zb — za for each term in the sum we recognize (4.8). It is then 
clear that this approximate solution can be extended to an exact solution in the vicinity of 
the wall, by correcting the locations Zi by a Taylor series in l/|zyi — zb\- 
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Eq.(4.7) shows that the solution to (4.8) exists as we approach the wall of marginal 
stability but ceases to exist as we cross the wall since the left hand side of (4.7) and hence 
\za ~ zb\ computed from (4.8) now becomes negative. Thus the jump in the index will be 
given by the contribution to (2.38) from this class of fixed points. To evaluate this jump, we 
first note that 

yT.i<j aij sign[zj-Zi] ^ yT.i.jeA,i<j sign[2j-2,] yT.i,jeB,i<j "ij sign[zj~Zi] y-fAB s\gn[zB-ZA] _ (4.11) 

Similarly, 

Finally we need to compute s{p) near the wall of marginal stability. This has been analyzed 
in detail in appendix A.l. The net result is that for the configuration of the type we are 
considering, s(p) is given by 

s{p) = sa sb sign(7yiij) sign[zij - za] ■ (4.13) 

where sa is the sign which appears in the contribution to ua centered collinear configurations 
with centers at {zi,i G A}, and is the sign which appears in the contribution to ub 
centered collinear configurations with centers at {zi,i G B}. 

Putting (4.11), (4.12) and (4.13) into (2.38) and summing over all n^i-centered solutions 
involving the charges in the set A and all n^-centered solutions involving charges in the set 
B we see that the net contribution from solutions of this type is given by 

-signi-fAB)sign[zB-ZA\ {y-y-y (-y)^-^^^^^"^!^^-^^! (?,ef({«.,^ e A};y)g^M(^„j G B};y) . 

(4.14) 

Adding up the contributions from the configurations for which Zb — za > and those for 
which zb — Za < 0, we get the net contribution to the index from the solutions which decay 
across the wall of marginal stability: 

-sign(7^B) {y - y'r' [i-yV-'^ - (-2/)"^^"] 9rA{(^^, ^ e A}; y) g,,,{{a,,3 e B}- y) . 

(4.15) 

This is the correct wall crossing formula for primitive decays. Besides the solutions considered 
here, (2.24) can have other solutions for which the relative distances between the centers 
remain finite as we approach the wall of marginal stability. They will continue to exist on 
the other side of the wall and must decay at other walls of marginal stability before we reach 
the attractor point. 

The analysis can be easily generalized to the case of general non-primitive decays. We 
shall sketch the derivation below. Suppose that we are near a wall of marginal stability where 
the total charge 7 can decay into L states carrying charges 7^ = XlieA «i for m = 1, . . . , L. 
For this we need the charges 7^ for different m to lie in a two dimensional plane and have 
^i&Am ~^ ^ each m. We approach the wall from the side where 

(7m, 7n) ( 5Z (4-16) 
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for any pair {Am, An). In this chamber there exists a class of solutions in which the elements 
of the set Am are bunched together for each m within a finite distance and the relative 
separation between the elements of the set Am and the elements of the set An go to infinity 
for every pair {m,n). These are the solutions which disappear across the wall of marginal 
stability; hence the change in the index is given by the index associated with these configu- 
rations. We shall order the sets Am such that (7m,7n) > for m < n. Let us assume that 
for a given collinear solution p of this type, the sets Am are arranged along the 2;-axis as 
Aa-(i), . . . , ^o-(L) for some permutation cr(l), . . . , cr(L) of 1, . . . , L. Then the contribution to 
the summand in (2.38) from such a configuration will be given by 

L 

S{p) Y[ (|/^«,jeA^"i.sig°fe-^«]) J-j- ^^(7a(™),7a(u)>^ _ (^417) 
m=l m<n 

The sum over collinear fixed points p which respect the bunching of the centers into the 
sets Ai, . . . ,Al will involve independent sum over collinear fixed points inside each set Am, 
generating the contribution grci{{i £ Am}', y) and the sum over permutations of the sets Am, 
generating a contribution grei{{lm}',y)- The final result, multiplied by the Hi ^ref ({'^i}' 
will be given by 



C/ref({7m};2/)n 



ieAr, 



(4.18) 



This is precisely what is needed to produce the wall crossing formula given in [17]. The 
gre{{{cii,i £ ^m}',y)Y[ieA^^Tcii{^i}'y) ^ctor Contribute to the index of the m-th 'black 
hole molecule' introduced in [17], while grcf{{lm}',y) is the function multiplying the indices 
of black hole molecules in the wall crossing formula of [17]. 

4.2 Consistency with split attractor flow conjecture 

The above analysis also shows that our proposal is consistent with the split attractor flow 
conjecture under certain assumptions. To see this let us consider the system of black holes 
carrying charges and let the moduli flow all the way to the attractor point. On 

the way we may cross several walls of marginal stability. If at the attractor point there 
are no multi-centered solutions, the result for the index at the original point in the moduli 
space can be computed by adding the contributions from the jumps across different walls. 
Suppose on the s-th wall the system is marginally unstable against decay into a pair of 
states carrying charges J2ieAs ^ieBs "^^^ Then the index associated with the states 

which decay across the wall is given by (4.15) with {A,B) replaced by {As,Bs). Thus the 
net contribution to the index at the original point from this multi-centered configuration is 
given by 

9reA{(^ihy) = - 2^ sign(7^^Bj _ 1 

s y y 

g^^^{{ai,ie As};y)g,,i{{aj,j e Bs};y) . (4.19) 
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In the y — 1 limit this gives 

9{{ai}) = $^(-ir^^"=+'l7A.Bj^7({«.,^ G A,})g{{a,,j G B,}) . (4.20) 

s 

We can now take the multi-centered configurations carrying charges {ai^i G As], {oij,i G 
Bs] and calculate their indices by fiowing along their attractor fiow lines in the same way. 
Continuing this process till we are left with only single centered black holes, we arrive at the 
split attractor fiow conjecture. 

Note however that the above analysis relies on an assumption: that the only jumps in 
g'j.gf (and in g) take place at the walls of marginal stability. In particular the solutions should 
not disappear away from the walls of marginal stability (or if they do then they disappear 
in pairs so that there is no net change in the index). 

5. Index from scaling mult i- centered solutions 

The goal of this section is to determine the contribution of scaling solutions to the index 
of a multi-centered black hole configuration using the Coulomb branch analysis. As shown 
in appendix B, in the absence of scaling solutions the functions G defined in (4.3) are 
automatically Laurent polynomials, as is required in order for the result to be a bona fide 
character of SU{2). This property however does not hold in general when scaling solutions 
are present, and (4.3) has to be corrected. We shall denote by G'^°f the contribution due 
to regular collinear fixed points only, and determine the corrections to (4.3) by requiring 
that after adding these corrections the final expression must be a proper SU{2) character, 
whenever the single-centered refined indices Q^{a,y) are SU{2) characters. 
Specifically, we propose to modify (4.3) into 

(^ref (7, y)= Yl G'eoll({A}, {m,}; y) J] (fif.f (A, y""^) + l^scaling(A, J/"')) , (5-1) 

{/Jier},{miGZ} i 

where riscaiing(«, y) is given by 

acaiing(a, y)= Yl ^(i/^^}' i'^*}' n ^-f 2/"') ' (5-2) 

{/3i6r},{™iez} i 

for some function if({/3j}, {m,j}; y) to be determined. To determine H we substitute (5.2) 
into (5.1) to express the latter equation as 

^reih,y)= E G({A},{m,};y) n^fcf(A,l/'"0, (5.3) 

{/3,6r},{m,eZ} i 

for some functions G. We fix H by requiring that G({/3j}, {mj}; y) are given by Laurent 
polynomials in y. This leaves open the possibility of adding Laurent polynomials to H. This 
ambiguity is resolved using the minimal modification hypothesis, which says that H must 
be symmetric under y — )■ y~^ and vanish as y — )■ oo. 
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In practice we solve for the functions H using an iterative scheme involving the number 
of centers. For this suppose we know {mj}; y) in all cases for 'Y^^mi < {n — 1). 

Now we can substitute (5.2) into eq.(5.1) and compute the coefficient of Ili^reflA'^™') ^^r 
all terms with Yli'^^i — The only unknown term is f^scaiing(7; y); originating from the 
replacement of fl^i'j, y) by ^'^(7, y) + f2scaiing(7, y) in the right hand side of (5.1). This gives 

i/({A},{m,};y) J] ^feflA, 2/"^') • (5-4) 

Thus requiring the coefficient of Hi ^rcf (A? y'"') to be a Laurent polynomial in y we can 
determine if({/3j}, {mj}; ?/) for ^iTrii = n. This procedure can then be repeated to find 
H{{Pi}, {'mi}; y) for rrii = n + 1 and so on. 

The algorithm given above gives a prescription for finding the net contribution from 
the scaling solutions for a fixed total charge 7. In the rest of the section we shall see how 
this prescription can be used to determine the contribution of the scaling solutions to a 
configuration containing a fixed set of centers. 

5.1 Correction to g^^i for non-identical centers 

We shall now show that if the centers carry non-identical charges then the minimal modifica- 
tion hypothesis translates to a simple rule for correcting the function (/ref (oi, • • • , 2/)- First 
suppose that ai, . . . , a„ have been chosen such that there is a scaling solution where all the 
centers come together, but no scaling solution where a subset of the centers come together. 
In this case our proposal for the contribution of the scaling solution to 5'ref(«i, • • • ,Cin',y) 
gives 

tog(«i,---,«n;i/) = (-l)^-^"-+"-'(y-y-')'-" Yl {y' - {-iry-'} , (5.5) 

0<fe<(n-2) 

where a^'s are constants to be adjusted so that, after adding (5.5) to the contribution (2.38) 
of the coUinear fixed points, the result has a finite limit as ?/ — )■ 1. It is easy to see that 
the number of a^s is precisely equal to the number of divergent terms in the y — )■ 1 limit, 
so that requiring finiteness as y — )■ 1 uniquely fixes all the a^s. For example for even n the 
possible values of k range from 1 to (n — 2), with k taking either only even or only odd values 
depending on the parity of YliKj'^ij- The number of a^s is then [n — 2)/2. On the other 
hand using y ■v^ y~^ symmetry we see that the divergent terms are of the form {y — y~^)~'^^ 
for s = 1, 2, . . . , (n — 2)/2, giving precisely (n — 2)/2 possible divergent terms. For n odd the 
allowed values of k are in the range to (n — 2) and again k takes either only even values or 
only odd values. This gives (n — l)/2 possible a^'s. On the other hand the possible divergent 
terms are of the form {y — y^^)^"^^ for 1 < s < {n — l)/2, giving precisely (n — l)/2 possible 
divergent terms. It is also worthwhile to note that the condition A; < n — 2 on the powers of 
y appearing in the numerator of (5.5) is equivalent to the requirement that the correction 
(5.5) vanishes in the limit y — )■ 00. 

Now consider the more general case where there are scaling solutions in which a subset 
of the centers come together. As mentioned earlier, we shall proceed by induction, i.e. 
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assume that all the fixed points (including scaling solutions) and their contributions have 
been determined for any number of centers less or equal to n — 1 and then show how this can 
be used to infer the result for n centered black hole solution. Let us consider an n-centered 
black hole configuration with centers carrying charges ai, . . . , a„. Now such configurations 
will include a set of 'regular' fixed points where all the centers are separated along the z-axis. 
We can determine them using numerical methods. To those we need to add the contribution 
from fixed points where a (subset of) the centers lie on top of each other. A generic fixed 
point of this type will have the charges {a, (o, a, (o, • • • , a^e)} lyi^ig ^op of ^ach other 
for 1 < / < s, > 1 such that Yli — ^- To determine its contribution we proceed as 
follows. We first consider an Z-centered configuration with individual centers carrying charges 
a (I) H a (I) with 1 < / < s and determine all its regular fixed points. In this computation 

the effective constant q for the Z-th charge will be given by 2Im(e~*"Za (jj+.-.o (;))• Now we 

multiply this contribution by the product of the weight factors of the s scaling solutions, 
with the rth solution containing centers {a, (o , a, (o , . . . , a, (o }. These are known by induction 
except for the case s = 1. We now add these to the contribution from the regular fixed points 
of the n-centered solution. This procedure leaves out the s = 1 term, corresponding to Z = 1, 
111 = n: this is the maximally scaling configuration, where all the centers are on top of each 
other. Our proposal for the contribution of this term to g^cf{c(i,- ■ ■ ,C(n]y) is again given 
by (5.5), where a^'s are constants to be adjusted so that after adding (5.5) to the other 
contributions (including the non-maximally scaling ones), the result can be expressed as a 
Laurent polynomial in y. 

The procedure just described generalizes the one introduced in §3.3 in the absence of 
scaling solutions. It relies on the assumption that the effect of the scaling solutions is 
to correct g^-cf into a proper SU (2) character, with the smallest possible powers of y and 
y~^ in the numerator (reflecting the classical fact that scaling solutions carry zero angular 
momentum). Eq.(5.5) together with (2.38) gives a complete prescription for computing the 
functions grci{{(^i}] y)- This construction guarantees that 5'rcf({ttj}; v) is given by a Laurent 
polynomial in y. As long as the charges are all different, the factor ^fref just described 
multiplied by YYi=i^^i<yi',y) determine the contribution to the index from multi-centered 
black hole solutions carrying charges ai,...,a„. The resulting expression is given by a 
Laurent polynomial in y. As discussed in §4, even when some of the Ofj's are equal but 
provided there are no scaling solutions, Eqs. (4.4), (4.5) give the correct result for the index 
of multi-centered black hole solutions. However when some of the a^'s are equal and scaling 
solutions are allowed, some additional corrections to (4.4) are necessary, which we shall now 
determine. 

5.2 Effect of identical particles 

We shall now consider the case where some of the centers carry identical charges, e.g. we 
have ri centers of charge (3i, r2 centers of charge /32 etc. We shall assume for simplicity that 
the only allowed scaling solutions involve all the centers coming together, - more general 
cases may be dealt with using the method of induction as before. In this case our proposal 
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for the index associated with this configuration is the following generalization of (4.4): 



Wez;nfe,4'')>i, (5-6) 



where the function /cor vanishes if all the r^'s are 1, but can be non-zero if some of the r^'s are 
larger than 1. The function g^-cf is defined as in (4.5) with g^ci including the corrections due 
to scaling solutions described in §5.1. The need for the additional correction terms /cor can 
be seen by noting that due to the presence of the {y^i — y"'*' ) factor in the denominator 
it is not guaranteed that the first term in (5.6) can be expressed as a Laurent polynomial in 
y. /cor is adjusted to compensate for this. We choose 

/cor({/34,{r4;2/) = E M{/3/^};{4"^};y)nn^'f(A,y^'''^ (5.7) 

{n.},{4"'} ' 
nk,s^k^ &^;'n-k,s^f°'^ >1, 

where h{{Pk}', {s^k^}] u) is chosen so that 

1. It is invariant under y — )■ y~^. 

2. iim,^ooM{/34;{4"^};2/) = 0. 

3. (5.6) has an expansion of the form Xlmez ^mU"^ with a finite number of terms whenever 
the f2^f(/3j, ?/)'s have this property. 

Note that the second condition above is another manifestation of the 'minimal modification 
hypothesis'. These three requirements fix the function h completely in any given situa- 
tion. Thus this prescription gives a complete algorithm for computing the spectrum of 
multi-centered black holes given the coUinear fixed point solutions to (2.24) satisfying the 
requirement (2.10). 

The term in which all the s^^^'s in the argument of h are 1 would combine with the 
term where each s^lf* = 1 in the first term in (5.6). Since the latter terms do not have any 
unwanted denominators and are automatically given by Laurent polynomials in there is 
no need for any correction terms. Thus /i(/3i, /32, ■ ■ ■ ; {-s^.'^''}; y) vanishes if all the s[."'''s are 1. 
This shows that if all the r^'s are 1, i.e. the a^'s are all different, then /cor vanishes. 

5.3 Illustration 

Since the above discussion has been somewhat abstract we shall now demonstrate this pro- 
cedure by a hypothetical example (which will map on to a real example in §6.5). Sup- 
pose we have a four centered solution with charges ai, Q!2, = ti/3, = r2/9 such that 
(ai, (3) = —{a2, (3). Suppose further that in some region of the moduli space we find that the 
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only coUinear configurations involve the order 1, 2, 3, 4 and its mirror 4, 3, 2, 1 with s{p) = ±1. 
Then the net contribution to (7rcf ■ ■ ■ ,C(A]y) from the collinear fixed points is 

(_l)"i2+3(^ - y-'r^y'"'' - . (5.8) 

Let us suppose further that there are no solutions where a proper subset of the centers 
are in the scaling configurations, - the only additional contribution comes from the scaling 
solution where all the centers approach each other. We can determine this term according 
to the prescription (5.5). This gives the net contribution to (^rcf (oi, • • • , «4; y) to be 

(?ref(«i, ...,a,;y) =(-l)°-+3(y - y-^y^ - - 1 a^2iy' - Z/~')} for ai2 G 2Z, 

= (_l)"i2+3(y - y-'y^ {y^^^ - y-"^^ - au{y - y'')} for ai2 G 2Z + 1 . 

(5.9) 

After taking out a factor of y^"^^ we see that the term inside { } is a polynomial in y^ and 
has triple zero at y = 1. Thus it must have a factor of (1 — y"^)^ cancelling the (1 — y^)^ 
factor in the denominator, and the expressions for 5'ref(tti, • • • ,tt4;y) given in (5.9) can be 
expressed as a Laurent polynomial in y. 

As long as all a^'s are different this ends the discussion for the contribution from this 
four centered terms. Suppose however ri = r2 = r so that as = a^. In this case the net 
contribution to the index from this four centered configuration will be given by 

^5(rof(ai, a2, rf3, r/3; y) fifef ("i! v) ^fcf ("2, y) ^rei{rf3, yf 

+ \ ^7rcf («i, «2, 2r/3; y) \ fif^^ (ai, y) fif.f (^2, y) ^L{rP, y^) (5-10) 

2 _ y-I 

+ /i(«i,a2,r/3;4'^ = 1,4'^ = 1,4'^ =2; y) fiL(ai,y) fifef(«2,|/) ni,{rl3,y^) . 

The second term comes from the term 5'ref(ai, ^2, 2r/3; y)f2^£(ai, ?/)f2^f(a2, y)VL^^,^{2r [3 , y) and 
the third term is the correction term given in (5.6), (5.7). To proceed we need to know 
g{ai, ^2, 2r/3; y). Suppose at the same point in the moduli space the collinear three centered 
configurations carrying charges ai, 02, 0^3 = 2r/3 are of the form 123 and 321. In this case 
we have 

(7rcf («i, «2, 2r/3; y) =(-l)°-+2(y - y~^y^ {y-- + y""- - 2} for G 2Z, 

= (_l)-i2+2^^ - y^Y^ + y-"- -{y + y^^ for ai2 G 2Z + 1 . 

(5.11) 

Note that we have added the correction terms due to the scaling solutions according to 
(5.5). We can now substitute (5.9) and (5.11) into (5.10) and determine h by requiring 
that the resulting expression is given by a Laurent polynomial in y for any choice of fi^f 
satisfying similar properties. Now the first term clearly has this property. In the second 
term (7ref (ai, 0^2, 2r/3; y) has this property, but the factor of {y — y~^) / {y"^ — y~'^) = y/il + y"^) 
has a factor of (1 + |/^) in the denominator which could potentially spoil this property unless 
the numerator has a factor of (1 + 1/^). In this case we need a non- vanishing h to cancel the 
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unwanted terms. To proceed we note that if ai2 is odd (ai2 = {2k + 1) with k E then we 
have 

I ^ grAc^i, «2, 2r/3; y) = -\y-''^\l + y')-\l - y'Y^l - y''){l - y''^') . (5.12) 

Each of the two factors (1 — y^^) and (1 — ^^^+2) has a factor of (1 — y^) cancehng the 
(1 — y^)^^ factor. Furthermore (1 — y^^) for even k and (1 — 7/^^+2) for odd k also has a factor 
of (1 + y"^) that cancels the (1 + y^)"^. Thus in this case (5.10) gives a Laurent polynomial 
in y without any h term and we can set h = 0. On the other hand for even (0^12 = 2/c 
with /c G Z) we have 

1 grAc^u «2, 2r/3; y) = ^ y-'^'+^l + y2)-i(l - y2)-2(l - y'^ ■ (5-13) 

2 i + 2 

(1 — y^^')^ has two factors of (1 — y"^) cancelling the (1 — y^)"^ factor. Furthermore for even 
k it also has two factors of (1 + y"^) killing the factor of (1 + y^)"^. Thus again in this case 
there is no need for any correction term and we can set h to 0. Finally for k odd we can 
express the right hand side of (5.13) as 

ly-'''-Hl + yr'{l + y' + ---y''-'V- (5.14) 

Now as y"^ — ?■ —1 this term approaches \ y{l + y"^)^^- Thus from (5.10) we see that the 
unwanted terms may be cancelled by choosing h to be negative of this term. This gives 

/i(ai,a2,r-/3;si^^ = 1,4^^ = 1,4^^ = 2;t/) = for e 2Z + 1 

= for ai2 G 4Z ^5 j^g-j 

= -\{y + y'^Y^ for «i2e4Z + 2. 

In §6.5 we shall realize this example in the case of dipole halo configurations. 
5.4 Wall crossing re-examined 

Given the modifications due to the scaling solutions, we need to re-examine the analysis 
in §4.1 on the compatibility of our prescription with the wall crossing formula. Rather 
than doing a detailed analysis, we note that our prescription for including contributions 
from scaling solutions affects the factors g^ciiWi.i G Am}]y) IliGAm ^ref('^i;?/) square 
bracket of (4.18), as each set Am may allow for scaling solutions. Thus, the index associated 
with individual black hole molecules will change. However, the factor grci{{lm}'iy) in front, 
which determines the jump of the index under wall-crossing, will not be modified since the 
elements belonging to different sets always remain separated. Thus, the contributions of 
scaling solutions do not affect the consistency with the wall-crossing formula. 

6. Dipole halo configurations 

In this section, we verify our prescriptions for a class of multi-centered configurations whose 
phase space and quantization is completely understood, namely dipole halo configurations [16, 
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29]. These consists of ra-centered configurations with two distinguished centers carrying 
charges ai, 02, with ai2 7^ 0, and n — 2 centers carrying mutually local charges Oa, such 
that aia = —oi2a^OLah = 0. Such configurations were analyzed in detail in [16, 29], whose 
notations we follow. A particular realization of this system, which we shall use to frame 
our discussion, is given by a D6 — DQ pair with n — 2 DO-branes orbiting around it. After 
reviewing the main features of such configurations in §6.1, we describe the 3-centered case 
in detail in §6.2 (an example of which was already analyzed in §3.3.1), and present a general 
proof of the minimal modification hypothesis in §6.3, based on recursion relations for the 
equivariant index. Similar recursion relations for the equivariant volume are presented in 
§6.4. Explicit results for 4 and 5 centers can be found in Appendix C. In §6.5 we analyze 
four centered configurations with two identical centers. 

6.1 Generalities 

We shall consider a system of n-centers, the first one of which represents a D6 brane with 
certain U(l) flux, the second one describes a L'6-brane with opposite U{1) flux and the third 
one onwards corresponds to DO-branes carrying q^, . . . ,qn units of DO-brane charges. Thus 
the first and the second centers carry opposite DQ and D2-brane charges but the same DA 
and DO-brane charges. In the language used in §3.3.1 we have: 

= {1, -U, -KUy2 - B/2A, -kU^Q - BU/24:), 
a2 = {-l,-U,KUy2 + B/24,-KUyQ- BU/24), (6.1) 
tta = (0,0,0,ga), 

leading to 

«i2 = / = 4Kt/^/3 + i3[//6, aia = -qa, a2a = qa, "af, = 0, 3 < o, 6 < n , 

(6.2) 

where we assume that /, qa are positive integers. The system described in (3.20), (3.21) is a 
special case of this with n = 3. We work in the chamber of the moduli space where 

Ci = /i, C2 = -/i - ?7 ^ qa, Ca = r]Qa for a = 3, ... , n, (6.3) 

a=3,...,n 

where /i is a positive constant and 77 is a small positive or negative number. As described at 
the end of §3.3.1, the = subspace describes a threshold stability wall on which a subset 
of the DO-branes can get infinitely separated from the rest of the system. But the index does 
not jump across this wall, since the symplectic product of the charge vector of the expelled 
DO-brane charges with the total charge vector of the D6 — DQ — DO system vanishes. 

On the Coulomb branch, the system is described by multi-centered configurations satis- 
fying the equilibrium conditions 

= V , + V — = -/^- 6.4 

ria r2a ri2 4^ ria 

a=i,...,n 

In particular, at 77 = the Z^O-branes lie either at infinity or on the plane equidistant to the 
DQ and DQ branes. As shown in [16, 29], for this type of dipole halo configuration the phase 
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space of solutions to (6.4) at r/ = is toric and given by a T" ^ bundle over the polytope 
{qa}) = {(m„m) : < m„ < , -j < m < j , j > 0} C M""^ , (6.5) 

where 

J = 1/2- J2 ^« (6-6) 

a=3,...,n 

is a linear function of the n — 2 variables m^. These variables parametrize the angle 6a 
between fia and fi2 via = qa cos 6a, while m parametrizes the angle ^ between fi2 and 
the z-axis via m = j cos 6*. Physically m represents the component of the angular momentum 
along z-axis. The coordinates {m,ma) together with coordinates {(p,4>a) £ [0,27r]'^~^ along 
the torus fiber provide a set of Darboux coordinates on A^^, 

tu = —dm A d0 — duia A d0a . (6.7) 

a=3,...,n 

Denoting (pa = <Pa — <y where ^a'i^a = and using (6.6), it is straightforward to check that 
(6.7) agrees with (2.14) with u = ga sin 6'a d6'a A dcpa- The equivariant volume (2.15) can 
be rewritten as 

^ciassicai({«i};?/) = (-1)^""+^ 1^ ^ dmg ■ ■ ■ dm„ / dme^'"", 

= (-1)^-"+^ / dm3---dm„^^"^^(^-'^"-3^"^H ,^1 

^ ^ ' 0<ma<qa in ^ , y 



(6.^ 



For brevity, we shall denote the r.h.s. of (6.8) also by 5(1; {ga}a=3,...,n; '^)- Moreover, we let 
S{I] {ga}a=3,...,n; z^) = wheuever / < 0. 

For ga < //2, the upper bound in '}2ma < //2 is never attained, and the phase space 
is compact. The equivariant volume can therefore be evaluated rigorously using localization 
with respect to J3. For > 0, the fixed points of J3 are coUinear solutions to (6.4), where 
a (possibly empty) subset A C {3, . . . , n} of the DO-branes lie on the segment between the 
DQ and DQ along the z axis, while the DO-branes in the complement B = {3, . . . , n}\A lie 
on the semi-infinite z-axis extending from the D6-brane to infinity. In the limit 77 — 0, the 
centers in A coalesce at the mid-point between the DQ branes, while the centers in B run 
off to infinity. The distance between the DQ and DQ is given by ri2 = (/ — 2 XlaeA ^a) I ^1 
which is positive for any subset A. The angular momentum carried by this configuration is 
J3 = (/ — 2^^g^ga) /2. In appendix A. 2 we show that the sign s(p) evaluates to (— 1)"-*, 
where ha is the cardinality of the subset A (the order among the centers inside the clusters 
A or S is irrelevant since they carry mutually local charges). Thus, the classical phase space 
integral (6.8) evaluates to 

5(/;{ga}a=3,...,n;i^) = (21ny)-"+l J](-l)"-+^-"+l (/-2E.eA5» + (_i)-iy-^+2E„,.5.) . 

(6.9) 
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This can of course also be obtained by direct evaluation of (6.8). Since the fixed points are 
isolated, according to the discussion in §2.5, the exact refined index is obtained by replacing 
2\ny hj y — y'^, leading to 

A 

(6.10) 

In contrast, when 1/2 < J2a1ay ^^e phase space is non-compact, as it has a scaling 
region represented by the boundary Yla'^^o. = 1/2 on which j vanishes, and all centers ap- 
proach each other at arbitrarily small distances. Such configurations are therefore invariant 
under 5*0(3). To compactify A^n, we include the boundary j = i.e. supplement the open 
polytope V with the lower-dimensional polytope 

n 

Q(/,{g4) = {(m„m) : 0<m, <g„, ^ = 7/2 , m = 0} . (6.11) 

o=3 

Denoting by Ai^^^^ the (2n — 4)-dimensional compact toric manifold built over the polytope 
Q, with torus fiber parametrized by 0^, we define the compactification of M.^ as 

Mn = MnyJMT\ (6.12) 

and require that 5*0(3) acts trivially on the boundary Ai^^^^. 

The equivariant volume of M.n can in principle be evaluated by localization. One class 
of fixed points is given by the same type of collinear configurations as described above (6.9), 
with the proviso that the subset A must be chosen such that ri2 > 0, i.e. 

E^'^<^- (6-13) 

a&A 

The contribution of these isolated fixed points takes the same form as in (6.9), with the 
restriction (6.13) enforced. We denote this contribution by 5'coii(/; {qa\a=?,,...,n] t^)- In addition 
to these isolated fixed points, already present in Ain, one also expects a contribution from 
the submanifold A4^^^^ inside Ain- Although this contribution can in principle be computed 
using a generalization of the Duistermaat-Heckman formula [38], we shall not attempt to 
compute it directly, and instead define 

'S'(/; {qa}a=3,...,n; I') = 5coll(/; {ga}a=3,...,n; '^) + A5(/; {qa}a=3,...,n; I') , (6.14) 

where A5(/; {q'a}a=3,...,n; t-') is the contribution of Ai^^'^. By a direct computation of the 
equivariant volume using recursion relations in §6.4, we shall be able to read off A5(/; {ga}a=3,...,n; J^)- 
Similarly, the exact refined index 5'ref({tti}; 2/), which we also denote by 5'(/; {qa}] J^), decom- 
poses as 

S{I; {qa}a=3,...,n; 1^) = ScoM] {qa}a=3,...,n; I') + A5'(/; {qa} a=3,...,n] ^) , (6.15) 

^^In this example, there are no scahng configurations where only a subset of the centers scale together, 
since the total charge carried by such a subset would be mutually local with respect to the remaining, 
non-scaling DO-brane centers. 
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where S'coii(/; {qa}a=3,...,n',^) denotes the r.h.s. of (6.10), with restriction (6.13) enforced on 
the subset A, and AS{r, {qa}a=3,...,n] v) denotes the contribution of the fixed submanifold 
M.^^'^^. As explained in §5, the minimal modification hypothesis determines the correction 
AS* uniquely by requiring that is a SU{2) character and AS" — j- as z/ — )■ oo. 

On the other hand, since the phase space of n-centered solutions is a toric Kahler man- 
ifold, it can also be quantized exactly. We refer to [16, 29] for the details of this procedure 
in the dipole halo case, and merely quote the result. When all the g^'s are distinct, so that 
the DO-branes are distinguishable, the exact index is given by 

(7ref({a.};y) = (-l)^--+^ J2 ^ (6.16) 

where [xj denotes the largest integer smaller or equal to x, and the sum over m G |Z is such 
that^^ m — |(J — n + 1 — S^^m^) G Z . As already mentioned above (6.15), we denote by 
S{I] {qa}] i^) the r.h.s. of (6.16), and let S{r, {qa}', i^) = whenever I < n — 1. Performing 
the geometric sum over m, we arrive at a sum of SU{2) characters 

S(/;W;.) = (-!/-"« E ^-M(I-n^2-2j:UmM 

^-^ smh V 

m3,...,m„ 

maeZ,0<ma<ga-l, 
Y.a'^a<\il-n+\)l2\ 

If some of the ga's coincide, and if all fi(aj)'s are set to one, the index is still given by (6.16), 
(6.17), with the additional restriction that the m^'s corresponding to identical particles must 
be distinct, and the expression must be divided by a symmetry factor k\ for every set of k 
identical g^'s. 

In §6.3, we shall give an inductive proof of the minimal modification prescription de- 
scribed below (6.15) in this class of dipole halo configurations, by establishing recursion 
relations for the exact index (6.17) and comparing them with the recursion relations obeyed 
by its regular part 5'coii(/; {^a}; ^)- Using similar methods, in §6.4 we shall also demonstrate 
a variant of this minimal modification prescription, which allows to recover the exact equiv- 
ariant index S{1\ {^a}; J^) from the knowledge of the equivariant volume S{1\ {qa)\ v). This 
prescription goes as follows: 

1. In the expression for Sil] {^a}; ^\ replace all z/'s which do not appear as arguments of 
hyperbolic functions, by sinhi^. Let us denote the result by Sil] {^a}; ^)- 

2. If S{1\ {qa\\ v) can be expressed as Laurent polynomial in y = ^ then this is the exact 
result for S{1\ {^a}; ^\ Otherwise we add terms which vanish as — ±oo to make the 
expression into a Laurent polynomial in y. This gives the exact result for S{1\ {qa\\ v)- 

Before going to the general proof of these statements, we shall illustrate them in the 
case of = 3 centers in §6.2. Explicit computations for n = 4 and 5 centers can be found in 

^^This condition was not stated explicitly in [16, 29] but is necessary for the consistency of the SU{2) 
action. 
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Appendix C. In §6.5 we shall study four centered configurations with two identical centers 
and use the minimal modification prescription of §5.2 to compute the index. This is then 
compared with the known exact results. 



6.2 Three-centered solutions 

For n = 3 we have to consider two cases separately. 

6.2.1 Non-scaling case: qs < ^ 

For ga < |, the sum over coUinear configurations (6.9) reduces to 

(-lY 



Scon{l; gs; y) 



(6.18) 



This can be expressed as a Laurent polynomial in y, and indeed agrees with the exact refined 
index computed from (6.16). 

Alternatively we can begin with the equivariant volume (6.8). This gives 

Sil;qs;y) = i-iy f dm, ^^"^^^^^ " ^"^'^'^ = (-i)^ ^"^^^^^) " " . (6.19) 

Jq V 1v 

After the replacement v sinh v in the denominator, we get 

rCOsh(Jz/) — cosh[(/ — 1qz)v\ 



S{I;qs;y) = {-iy 



2 sinh u 



(6.20) 



It is easy to see that this agrees with (6.18) and hence can be expressed as a Laurent 
polynomial in y. Thus in this case there is no need to add any correction terms and S'(/; q,; y) 
is the same as 5(1; gs;?/). In particular, the symplectic volume 5(7; gs;?/ = 1) agrees with 
the exact index S'(/; ga; y = 1) in the non-scaling regime. 

6.2.2 Scaling case: gs > | 

In this case (6.9), with the restriction (6.13) on the set A, reduces to 



Sco\\{l; q3;y) 



(y-y^^y 



y' + y ' 



(6.21) 



This diverges as ?/ — ?■ 1 and hence we must add corrections described in (5.5). This indeed 
reproduces the exact index computed from (6.16), 



S{i;q3;y) 



{y-y-'y 
{y-y-'y 



y^ + y ^ - (y + y ^ 
y' + 2/"' - 2 



for J odd 
for / even . 



(6.22) 



On the other hand the result of the equivariant volume is insensitive to the parity of /, 
and gives 



5(/;g3;y) = (-!)' 



1/2 



dm' 



sinh[(/ — 2m,)v\ 



■1) 



jCosh(/z/) — 1 



21^2 



(6.23) 
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The last term in the numerator is recognized as the contribution of the scahng fixed point 
Q = {rris = I/2,m = 0}. Replacing u — )■ sinh u in the denominator we arrive at 



, j Cosh(/t^) - 1 



2 sinh V 



(6.24) 



For / even, this agrees with the exact result (6.22). On the other hand if / odd, S differs 
from the exact result S by 



J cosh V 



2 sinh^ u 



for / odd . 



(6.25) 



Note that the correction term vanishes as z/ — ^ ±oo. Thus, had we had started with the 
expression for S and added corrections following the prescription given below (6.17), we 
would have arrived at the correct expression for S. It is also worth noting that when I is 
odd, the exact index S'(/; gs; y = 1) differs from the symplectic volume S{I] q^, y = 1) by a 
fraction (—1)^/4, which is necessary to make the result integer. Thus the non-renormalization 
property which was observed in the non-scaling case breaks down in this case. Finally, we 
note that we have not discussed the effect of the regularity conditions (3.8). This was 
discussed in a special case in §3.3.1, and we expect this to be satisfied also for all the 
solutions described in this section. 



6.3 An inductive proof of the minimal modification hypothesis 

In this section, we shall prove by induction that AS {I; {qa}a=3,...,n', J^) appearing in (6.15), giv- 
ing the difference between the exact index (6.16) and the contribution Scon{I', {qa}a=3,...,n] ^) 
of regular collinear configurations given in (6.10), vanishes as z/ — )■ oo. This proves the 
validity of the minimal modification hypothesis for all n for the specific system under con- 
sideration. 

First, we establish recursion relations for the equivariant index S*, defined in (6.17). For 
this purpose it is useful to introduce a variant of S defined by 



C{I] {qa]a=3,...,n\ 1^) 



\I+n-l 



E 



cosh 



(/_^ + 2-2E,=3,...,n"^a) 



V 



m3,...,m„ 

ma&,0<ma<qa — l 
T.a"^c.<W-n+l)/2\ 



sinh v 



and another quantity 



E{I-, {qa] 



a=3,...n ) 



J+n-l 



E 



1 . 



(6.26) 
(6.27) 



m3,...,m„ 

ma£l',0<ma<qa-l, 

E.™-=La-"+i)/2J 



This is recognized as the number of integer points closest to the polytope Q in (6.11). We 
also define 5", C and E to be zero for I < (n — 1). For n = 2, i.e. in the absence of DO-branes, 
we let 

(6.28) 



sinh u 



sinh u 
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if J > 0, and zero otherwise. The quantity E{I] {qa}) can be evaluated from (6.26) using 

{qa}a=3,...n 



hm 



V {C{I - 1 - 2g„; {qa}a=3,...,n-i; l^) - C{I - 1; {qa}a=3,...,n-~i; T^)) 

(6.29) 

vahd for 72 > 3. 

Performing the sum (6.17) over the last charge m„ first, we see that the sum is empty if 



X]a=3 "^a > [(/ — n + 1)/2J , or else runs from to g — 1 where q = min(q'„, [(J — n + l)/2j 



X]a=3^a + -'-)• Using the geometric sum identities 

q-l 

■ -i / A o n N cosni/i -h \zq 

2 sinh B 

smh[A + {2q-l)B]- sinh(v4 - B) 

josni^/i -t- An ni) = 

m=Q 



• t,/. nr. ^ cosh[A+{2q-l)B]-cosh{A- B) 
> smh(A + 2Bm) = — ^—^ ^ 

m=0 
9-1 

cosh(A + 2Bm\ 



2 sinh B 



(6.30) 
(6.31) 



with y4=(/ — n + 2 — 2 X]a=3 "^a)^ ^^"^ B = —u, we find, for any n > 3, 

1 



'S'(/; {ga}a= 



C{r, {qa}a- 



■3,...,ni 



V] 



2 sinh V . 
cosh^l + (-l)^-"sinh2 ^ 



:3,...,ni 



2 sinh^ u 
1 

2 sinh 1/ - 

1 H 



C(/ - 2g„; {ga}a=3,...,n-i; Z^) - (^(Z; {qa}a 
v_ 

~^E[1- {g4a=3,...,n) 



S{I - 2g„; {ga}a=3,...,n-i; J^) " 'S'(/; {ga}a=3,...,n-i; Z^) 
iM-n 



(6.32) 



4 sinh z/ 



-E{I; {qa}a- 



-3,...,n 



Similarly, we can derive a recursion relation for the contributions of the regular collinear 
solutions to the refined index. Given a collinear configuration contributing to the (n — 1) 
particle system with the DO-branes carrying charges q^, . . . , qn-i, we can construct a collinear 
configuration contributing to the n particle system with an additional DO-brane with charge 
qn as follows. For simplicity we shall work in the ?7 — )■ limit. First of all the additional 
Z^O-brane can always be placed at infinity for any collinear configuration of the original 
{n — 1) particle system. Also, for any collinear configuration of the original system, if 
the set A containing the DO-branes at the midpoint between D6 and D6 brane satisfies 
XlaeA la < {I ~ '^qn)/'^i wc cau add the ra-th DO-brane at the midpoint. Taking into account 
the various signs appearing in (6.10) we get the recursion relations, for any n > 3, 



'S'coll(-^; {qa]a=3,...,n-, Z^) 



Ccoll(-^; {qa]a= 



-3,...,ni 



2 sinh u 
1 

2 sinh V 



Cco\\{I — '^q-n', {qa}a= 



:3, 



-i; u) — Ccoll(/; {qa}a=3,...,n-l', J^) 



'S'colK-^ — 2g„; {qa}a=3,...,n~l', ^) — Sco\\{I\ {qa} a=3,...,n-~l'i Z^) 

(6.33) 



where C*coii is defined with an opposite sign compared to S'coii in (6.10), 

C'coll(/;{ga}a=3,...,n;Z^) = (^-^ ')""+' ^("1)"""''^""' (/^^E^e... _ (_ E.e. -Z^) 

A 

(6.34) 
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Subtracting (6.33) from (6.32), we arrive at 

AS{r, {qa}a=3,...,n; l^) = ^ . AC'(J - 2g„; {qa}a=3,...,n-l] " AC'(/; {qa] a='i,...,n-l\ 



2 sinh u . 
cosh^^^ + r-n^-^sinh^!^ 



o • ,2 -E{I; {ga}a=3,...,nj 

zsmii z/ 

AC'(/; {ga}a=3,...,n; Z^) = o ■ , AS'(/ - 2g„; {ga}a=3,...,n-i; T^) " A5(/; {ga}a=3,...,n-i; I^) 



2 sinh . 

1 + (-1)^-" 



E{I\ {ga}a=3,...,n) 



4 sinh z/ 

(6.35) 

Assuming that AS*, AC* -t- as ?/ — )■ oo for n — 1 particles, it immediately follows from (6.35) 
that the same statement holds for n particles. The validity of the assumption at n = 3 
is easily checked using the explicit results in the previous subsection. Thus, the minimal 
modification prescription is proved for this class of multi-centered configurations. It would 
be interesting to compute A^* directly by using the formula (2.35). 

6.4 Recursion relations for the equivariant volume 

In this section we shall derive recursion relations for the equivariant volume (6.8) similar to 
the ones given in §6.3 , and then use them to prove the prescription given below (6.17). We 
first define a variant of the equivariant volume (6.8), 

C{I; {m -) = (-1)^+"-^ dm3 ■ ■ ■ dm„, eosh[(/ - 2 mj.] ^^^^^^ 

and 

» n 

E{I-{qa]a=3,...,n) = (-1)'+"-' / dmg " " " dm„ 5 ( ^ - //2) . (6.37) 

J 0<ma<qa „_3 

This last expression is recognized as the symplectic volume of the submanifold of fixed points 
based over the polytope Q. For = 2, i.e. in the absence of DO-branes, we set 

s{i. V) = (-1)^-^ ""^"^ , ai, u) = (_i)/-i ^"^M/.) 

if / > 0, and zero otherwise. E can be evaluated in terms of (6.36) using 

E{I\ {ga}a=3,...,n) = lim [v {C {I - 2g„; {qa}a=3,...,n~i; - C{I] {qa} a=3,...,n-i; ^))] , (6.39) 
1/— >0 

or directly from its Fourier representation, 

E{I;{qa}a=3,...,n) = (-1)^+""^ / ^T^e-^^^-S-^-,-." J] 2 siu ft . (6.40) 
In particular, it is a piecewise polynomial in /, g^. 
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By integrating first with respect to the last charge rrin in (6.8), (6.36), it is straightfor- 
ward to estabhsh the following relations, valid for n > 3: 

S{I; {qa}a=3,...,n; V) =^ {C{I - 2g„; {ga}a=3,...,n-i; v) - C{1\ {qo) a=Z,...,n-\\ v)\ 

^ (6.41) 

- ^ {ga}a=3,...,n.) 
C(^; {ga}a=3,...,n; ^) = ^ ~ ^^n! {Qa} a=3,...,n~i; l^) - S{I] {qa} a='i,...,n-l] T^)] (6.42) 

Similarly as in (6.33), the contribution to S{I] {qa}) from the regular, coUinear fixed points 
satisfies 

'S'coll(-^; {qa\a=Z,...,n'-, T^) [C'coll(-^ — '^q-n'-, {9a}a=3,...,n-i; ^) ^ C'coll(-^; {Q'a}a=3,...,n-i; T^)] 

Cco\\{I'-, {qa}a='i,...,n'-, l^) [Sco\\{I — '^qn] {qa} a=3,...,n-l', i^) — 'S'coll(/; {qa} a=3,...,n-l', J^)] 

(6.43) 

where Ccou is defined with an opposite sign compared to S'cou in (6.9). Applying (6. 41), (6. 43) 
recursively, the difference AS = S — S'cou evaluates to 

A^(/; {qa}a=3,...,n, ^) = 2 ^ E L)n-"^^^ ^^^ " 2 ^ {qa} a=3,...,.^) (6.44) 

ri — m£2Z 

where for given m, the sum runs over subsets B C {m + 1, . . . n} subject to the restriction 
that Ylih<^B'ii' ^ -^/2- This difference can be understood as the contribution of the fixed 
submanifold J^^^^^ built over the polytope Q. In particular, the first term with m = n in 
the sum (6.44) is just the symplectic volume E{I, {qa}a=3,...,n) of J^n^'^, rescaled by a factor 
— l/(2i^^). The remaining terms in (6.44) should originate from the Euler class of the normal 
bundle to Ai^^'"^^ which determines the integration mesure on the fixed submanifold, as in 
(2.35). We give evidence for this claim in two specific examples with = 4 and = 5 in 
§C.3. 

We now turn to the proof of the prescription given below (6.17) for obtaining the exact 
equivariant index from the equivariant volume. Recall that S was defined as the result of 
replacing in S all the i^'s outside the argument of hyperbolic functions by sinhz/. We define 
C as the result of a similar replacement in C. Then the recursion relations for S and C can 
be obtained by modifying (6.41), (6.42) to: 

S{I] {qa}a=3,...,n] l^) = ^ .\ C{I " 2g„; {qa}a=3,...,n-l] 1^) " C{I] {qa} a=3,...,n-l] l^) 



2 sinh V 
1 



- . 2 -^(-^' {qa}a=3,...,n) 

1 smh V 

{ga}a=3,...,n; = ^ ^.^^ ^ S{I - 2g„; {ga}a=3,...,n-l ; " S{I\ {qa} a=3,...,n-\\ v) 



(6.45) 
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Subtracting (6.33) from (6.45), and defining AS* = S — Scou, AC = C — Ccoii, we arrive at 

1 



AS'(/; {qa}a=3,...,n; T^) 



2 sinh V . 
1 



AC(/ - 2g„; {ga}a=3,...,n-i; - ^C{I; {qa}a=3,...,n-i; ^) 



r, . 1 2 -^(^i {?a}a=3,...,n) 

2 smn u 



^C{I- {qa}a- 



:3,...,n) 



2 sinh 1/ 



AS'(/ - 2qn] {ga}a=3,...,n-i; «^) " AS'(/; {qa] a=3,...,n-l] I') 

(6.46) 



Assuming that AS*, AC — )■ as — oo for n — 1 particles, it immediately follows from 
(6.46) that the same statement holds for n particles. The validity of the assumption is easily 
checked for = 3. On the other hand we have already seen in §6.3 that 5* — Scow and 
C — Ccou vanishes as — ?■ oo. This shows that S — S and C — C vanishes as z/ — )■ oo, thereby 
confirming the prescription described below (6.17). 

6.5 Four-centered case with two identical centers 

So far we have only considered the cases where the centers carry distinct charges. For the 
dipole halo configuration analyzed here, the first case of identical charges arise for four cen- 
tered configurations in which q^ and q^ coincide. We shall now examine this case and compare 
the results obtained using the minimal modification hypothesis with the exact results. In 
this case our formula (5.6), (5.7) for the total index associated with this configuration is: 

1 



-Qreiiau "2, "3, ^3; y)fircf ("l. V) ^ref ("2, V) ^rcf ("3, VY 



1 

+ 7Tfi'rcf(ai,"2,2a3;?/) 



y-y' 



^fcf("i> y) ^fcf ("2, y) f^fcf ("3, y^ 



+ h{ai,a2,rP;s^^^ = 1,8^^^ = l,sf^ = 2; ?/)fif^f (ai, y) Q^^f{a2,y) (^f^f (as, y^) . 



(6.47) 



We shall determine the function h using the prescription of §5.2. Following the convention 
of [33] and appendix C we shall consider three separate cases A, B and D. The case C, for 
which ^3 < //2 < g4 is not relevant here since we have q^, = q^- 

Case A 

In this case we have 2^3 < J/2. Thus according to (6.18) we have 

-1)^ \ r r r . , . „ J (- 1 ) ^ ( 1 - ^ ^'^3 ) _ ) 



5'ref(ai,a2,2a3;y) 



{y - y^^Y 



y'-y'-'"' 



{y-y 



(6.48) 

Since the first factor has a factor of (1 + y"^) our prescription of §5.2 tells us that we have 
= 0. Thus using (6.47), (C.l) and (6.48) we see that the full index is given by 



2{y-y~^Y 
1 f 



-/+4g3 



^fef ^fef («2, y) l^fef ("3, y)' 



+ 



2{y-y~^Yy^-y^^) 



y -y 



-/+4g3 



ni,{ai, y) l^fef («2, y) J^fef («3, 1/') • 



(6.49) 
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It is straightforward to check that Eq. (6.49) with = 1 agrees with the exact 

refined index computed from (6.16), with the additional restriction rris < on the sum to 
account for the identity of the particles. 



Case B 

In this case we have < 1/2, 2q^ > 1/2, and hence according to (6.22) we have 



gvef{ai,a2,2a3]y) 



{y-y ^y 
{y-y~^Y 



y' + y ' -{y + y ') 



for / odd 



y' + y-'-2 



(6.50) 



for / even . 



This is the same as in (5.11) with ai2 replaced by /. Since h is determined from (^ref (^i, ^2, 2a3; y), 
it will be given by (5.15): 



h{ai,a2,rj3]s^l^ 



1 s^^^ 



1 S^^^ 



2; y) = for / G 2Z + 1 
= for / G 4Z 

1 ' 



{-iy+\y + y-'Y' for / G 4Z + 2 . 



(6.51) 



Using (6.47), (C.5), (6.50) and (6.51) we now get the total contribution to the index from 
this configuration to be 



1 (-1) 



/-I 



2{y-y-^f 

^fcf («b y) ^fef («2, y) ^ii{a3,yy 

1 {-ly 



+ 



2{y~y^^){y'-y-^) 



y' + y ' -{y + y ') 



fliiiai, y) fifef ("2, y) f^fcf ("3, y'^ 



for / G 2Z + 1 
(6.52) 



1 



'1) 



i-i 



y' - 2^^-^"^ 



y 



-/+2(73 



2 [y-y^^Y 

^fef ("1. y) ^fef ("2, y) ^fcf ("3, yf 

1 (-1)^ 



;(y'-r')(4g3-/) 



+ 



2[y-y-^W~y^^) 



y +y 



-2 



^ref ("!> Z/) ^ref («2, y) ^ref («3, 1/ 



for J G 4Z 

(6.53) 
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1 



-1) 



/-I 



2 {y - y-^f 
1 



-'^''■'-\{y'-y-'){^q,-i) 



+ 



^fef ("1, Z/) ^fef(a2, fifcf ("3, Z/)' 



2 



y' + y-' -2 



^ref(«l>l/)^ref(«2,l/)f^ref(«3,l/ ) 



for / e 4Z + 2 . 

(6.54) 

Again, one may check that Eq. (6.52)-(6.54) with fifef("'l/) = 1 

agree with the exact refined 

index computed from (6.16) with the restriction ms < m/^. 
Case D 

In this case we have ^3 > 1/2 and as a resuh g-[ci{(^i^oi2i2a^^y) and hence h will be given 
by (6.50) and (6.51) respectively. Thus the only difference from case B is in the expression 
for (jfrcf (tti, «2, ^3, as; y) given in (C.14). Thus using (6.47), (C.14), (6.50) and (6.51) we now 
get the total contribution to the index from this configuration to be 



-1) 



I -I 



2 {y-y-'f 
1 (- 



+ 



-y ^-Hy-y ^) 

1)^ 



fifrf(ai, y)nii{a2, 2/)l^fef («3, y? 



2{y~y^^){y'-y~^) 



y' + y~' - (y + y~^] 



for / G 2Z + 1 
(6.55) 



1 



-1) 



7-1 



2 (y-y-^r 
1 ( 



y^ -y ^ - ^{y^ -y ^) 



2 {y-y-'W-y-') 



y' + y~'-2 



^fcf ("1. 2/)^rrf(a2, 2/)fifef ("3, yf 



for J G 4Z 

(6.56) 



1 (-1) 



2 (y - y-^r 
1 ( 



/ - - ^ iy' - y~') 



+ 



2 (y - y-'){y' - y-') 



y' + y~'-2 



^fcf 2/)^fcf ("2, y)l^fef("3, y^ 

^fef ("1. 2/)^fef(a2, 2/)fifef ("3, y^) 



+ ^ (-l)^+=^(y + i/-^)-^fif,f (ai, 2/)fifef(«2, 2/)^^?ef («3, 2/') , 



for / G 4Z + 2 . 

(6.57) 



Again, one may check that Eq. (6.55)-(6.57) with Q^^^{a, y) = I agree with the exact refined 
index computed from (6.16) with the restriction < 7714. 
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A. Sign rules for collinear fixed points 

In this appendix, we provide some details on the computation of the sign of the contribu- 
tions from collinear fixed points near a wall of marginal stability (§4.1) and for dipole halo 
configurations (§6.1). Recall that s{p) is given by (2.29), where M[p) is the Hessian of the 
'superpotentiaF (2.25) at the critical point p. 

A.l Sign rules near a wall of marginal stability 

First we shall compute the sign s{p) associated with a collinear configuration p near a wall of 
marginal stability, if the collinear configuration breaks up into two widely separated clusters 
as we approach the wall. Such a configuration has been described in §4.1 where the two 
sets into which the centers split have been called A and B. We follow the notation of §4.1, 
and work with the configuration for which zb > za, i.e. zj — Zi > for i E A, j E B. Near 
R=\za — zb\ — oo, the superpotential (2.25) decomposes as 

W{zi, A) ~ Wa [zi^A, = — a) + Wb (^Zj^B, Ab = — a) - aij In \zi - Zj\ . (A.l) 

For such a configuration the Hessian of W with respect to A, ZifzA, -Zjes takes the form 

M = Mo + Ml , (A.2) 
where Mq is the Hessian in the strict i? — > oo limit and Mi is of order We have 

Mo = I ^UA Ma I (A.3) 

^ub Mb 



n 



where Ma is the Hessian of Wa with respect to Zi^A, ua is the n^-dimensional column vector 
with entries l/n^, and similarly for M^, ub- In particular, we note that Ma ua = MbUb = 0. 
Ml is given by 

(Mi)oo = (Mi)o* = (Mi)io = for 1 < 2 < n, 
''^ii Efcei?«ifc ior ij E A 

^ij Y^k^A^ki for i,j E B (A.4) 
—aij for i E A, j E B 
—Oiji for i E B, j E A 



(Ml).. 



1 
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On the other hand, the Hessian of Wa (respectively, Wb) with respect to A^, Za (a G A) 
(respectively, Xb, Zf,, b E B) is given by 

M. = (°r;). Ms = (°;;). (A.5) 

\ua MaJ \ub Mb J 

To compare the signs of det M and det Ma det Mb , we shall construct an eigensystem of 
M in terms of the eigensystems of Ma and Mb, in the limit i? — )■ oo. First, we note that an 
eigensystem of Ma is given by the + 1 (eigenvectors, eigenvalues) 

--^,ua],±^) , ([0,t;«],A«) (A.6) 

where f^'*, z = 1,...,?t,a — 1 are eigenvectors of Ma in the subspace orthogonal to the 
null eigenvector ua- Similarly, let v^^ be a system of eigenvectors of Mb in the subspace 
orthogonal to the null eigenvector ub, with eigenvalues A^'*. In the strict R ^ oo limit, M 
reduces to Mq and an eigensystem of Mq is given by the n^+riB + l (eigenvectors, eigenvalues) 



''^ua,'^ub],±^), ([0,v'-2\0b],X^^) , ([0,OA,t;g^],Ag')), ([0, n^, -mb], o) 

(A.7) 

In particular, the last eigenvector, corresponding to a change in the relative separation 
between the two clusters keeping their inner structure fixed, yields a zero- mode of Mq. Since 
the eigenvalues A^^^ A^^ are generically distinct and non-zero, the structure of the spectrum 
will retain its form at large but finite R, except that the last eigenvalue will be lifted to a non- 
zero eigenvalue X{R). To compute X{R), it suffices to use non-degenerate perturbation theory 
and calculate the expectation value of Mi on the unperturbed eigenvector U = [0,^^, —ub]- 
This gives 

^ ' UT-U uaUbR^ ■ ^ ' 

As a result, 

det M ~ -^A(i?) H A?) H Ag^ ~ --^ 7ab det Ma det Mb (A.9) 

« j 

Using (2.29), we arrive at (4.13) for zb > za- The result for zb < za follows by exchanging 
A and B in the above analysis. 

A. 2 Sign rules for collinear dipole halos 

In this section, we establish the sign rule s{p) = (—1)"-'* used in (6.9) for collinear dipole 
halo configurations in the non-scaling regime. For small positive t], the solutions described 
above (6.9) satisfy 

^ lizi + Z2) - VaGA, Zb lizi + Z2) - R \/beB, (A.IO) 



n n 
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where R ^ \/rv2hl- In this hmit, the Hessian M takes the following form to order R 

1 . . / 

Moo = 0, Mo, = Mio = - for 1 < j < n, M12 = M21 = --. — ^ , 



3. 



n [ri2y 



fees aeA "~ fees 



4 2 1 ri2 

Mia = Mai = J^r^^lcL + -^Qa for fl G A, Mifo = Mfoi = --^^ - -^Qh for 6 G 5 , 

4 2 1 ri2 

M2a = Ma2 = J^T^Ia " for a G v4, M26 = Mf,2 = -^Qb " -^Qh for 6 G 5 , 

Mafe = Mfea = for a G A, 6 G 5 . 

(A.ll) 

To evaluate the determinant of this matrix in the large R limit we proceed as follows. First 
we define a new matrix M by dropping the first two rows and columns of M; we have 
seen earlier that det M = — det M = — det M. To evaluate det M, we can add half of the 
the second to last row of the matrix to the first row and then add half of the second to 
last columns of the matrix to the first column. This does not change the determinant but 
simplifies the matrix. Let us denote the resulting matrix by Mo + Mi where Mq is the limit 
of the matrix as i? — )■ oo and Mi is the remainder. It is straightforward to check that Mo is 
diagonal and has eigenvalues 

a-2EaeA9a) / 8qa 



12 



{On,}, (A.12) 



where 0^^ denotes that the eigenvalue is repeated ub times. When we take into account the 
effect of Ml, the non-zero eigenvalues are not affected appreciably but the zero eigenvalues 
are lifted and can be obtained using first order degenerate perturbation theory. This gives 
the approximate eigenvalues of Mo + Mi to be: 



a-2EaeA?a) r 8qa \ f 2r 



12 



12 



ri2V r \\R' 



(A.13) 



To leading order in the limit i? — )• oo, det M = — det M is therefore given by 

det M - 2"+^"-^'^ ( (-^ ~ 2 ZlggA gg) (HagA gg) illb&B lb) fA 14) 

^ ' _^3(n-2-nA)^3nA-n+4 " \ ■ ) 

In particular, the sign of det M is (—1)"-'*+^. Note that to obtain this result, it is important 
to keep all subleading terms through order 1/ R^ as indicated in (A.ll), since the determinant 
at lower order vanishes. Using (2.29) we now arrive at 

s[p) = [-lYA. (A.15) 
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B. Laurent polynomial property in absence of scaling solutions 



Eq.(4.1), (4.2) (or equivalently (4.3)) gives the index associated with a multi-centered black 
hole solution when the charges of the components do not allow for scaling solutions. Let 
us restrict to the case where 7 is primitive so that (4.1) directly gives the refined index 
^rGf(7,l/) rather than its rational counterpart r2ref(7,i/). In this case the right hand side of 
(4.1) must be a Laurent polynomial in y (i.e. a finite linear combination of y^^) whenever 
the f2^f(aj, ?/)'s are since otherwise the result cannot be interpreted as an SU{2) character. 
Our goal in this appendix will be to prove this property of (4.1).^^ 

Clearly this will not be true for an arbitrary choice of the functions g^^^, but we shall 
use (4.19) - valid when at the attractor point only single centered black holes contribute 
to the index - to restrict the form of g^ci. It is clear from this equation that the right 
hand side of (4.19) will be a Laurent polynomial in y if the g^^^^s appearing on the right 
hand side have this property. Thus by iterative application of (4.19) we can establish that 
(/ref(tti, • • • , ttn, u) is givcu by a Laurent polynomial in y. It now follows from (4.1) that when 
the ttj's are all primitive then the right hand side of this equation is a Laurent polynomial in 
y since the individual Of^f (a^, y) = fif^f (a^, y)'s and the (^ref (^i, • • • , ctn, u) have this property. 
This argument fails when some of the a^'s are not primitive: in this case VL^^^{ai, y) defined 
in (4.2) have extra factors of mly"^ — y~^) in the denominator, which must cancel in order 
that the final expression is a Laurent polynomial in y. We shall now demonstrate that this 
cancellation does take place. 

By iterative application of (4.19) we can express 5'ref(«i; ■ ■ ■ 1 «n) as a sum over attractor 
flow trees in which a total charge 7 decays via an appropriate tree to the charges ai, . . . , 
At a vertex at which a charge /3i decays into two clusters of charge [32 and /^s, we get a 
multiplicative factor of 

(-l)<''-'««.ign(ft.A)^^^^ti^%^. (B.l) 

smh V 

Let us focus on the vertex at which a non-primitive charge ai gets attached to the tree. 
Suppose at this vertex an internal line carrying charge P + Ui decays into another internal line 
of charge (5 and the external line of charge a^. Our goal will be to show that the contribution 
from the vertex factor cancels the {y"^ — y~"^) factors in the denominator appearing in (4.2). 
The product of the vertex factor and the ri^{ai, y) factor is given by 

(-l)<'->«sig„(,3.a,)?''^'?«^n-(a..y) 

smh u 

= (-l)</^-)+^sign(/3, a.) Y: m-''^^^^^^^^Q^{a,/m, y-) . ^^'^^ 

^-^ smhfmi/j 

m\on 

Since m|Q;j, {fi-,o.i) is an integral multiple of m. In this case the unwanted denominator 
factors cancel and sinh((/3, a^)//)/ sinh(mi/) is a Laurent polynomial in y. 

^•^In order to prove that the result is indeed an SU{2) character, one must also show that the coefBcients 
of j/^"* are integers. This is indeed true, but we shall omit the proof. 
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This leaves the end vertices of the attractor flow tree, at which a charge + aj decays 
into a pair of charges a, and aj. The contribution from such a vertex will be of the form 



.l)(--.)+isign(a„ a,) 

smh 



■l)<"""^^+^sign(ai,aj) 2^2^m"V^ v\ j/ ; ns^^, „,m^ns^ 



sinh(mi/) sinh(pz/) 



^]^(«i/m,^/™)^]^(a,/p,l/^'). 

(B.3) 



Let us define z = y"^ = . Since m|aj and p|aj, {ai,aj) is an integral multiple of mp. In 
this case sinh((aj, ctj)//) will have zeroes at 2"^^ = 1. On the other hand the denominator 
smh{mv) sinh(pz/) has zeroes at z"* = 1 and also at = 1. If m and p are relatively prime 
then the locations of these zeroes are distinct except for a common zero of both factors at 
z = 1. Furthermore they coincide with the zeroes of the numerator at = 1. Thus all 
the zeroes of the denominator sinh(mz/) sinh(pz/) cancel against the zeroes of the numerator 
sinh((aj, sinhi/. Thus as long as m and p are relatively prime, (B.3) is a Laurent 

polynomial in y. 

Now suppose that q = gcd{m,p) > 1. Then sinh(mz/) sinh(pz/) will have double zeroes 
at each of the q solutions to z"^ = 1. In contrast the sinh((aj, factor in the numerator 

will generically have only single zeroes at each solution of 2;'' = 1. Combining this with the 
extra factor of sinhz/ in the numerator, we see that there is effectively a left-over factor of 
sinhi// sinh(gi/) from the vertex, multiplied by factors which are Laurent polynomials in y. 
We now show that the factor sinh u/ sinh(gz/) is cancelled by other vertex factors in the same 
tree. 

To see this, note that since q divides both ctj, aj, it divides their sum + aj. We can 
now repeat the analysis at the next vertex where say a line carrying charge at + aj + (3 
splits into charges /3 and + aj. As long as (3 and q do not have a common factor, the 
analysis of the previous paragraph shows that the vertex factor will cancel all the unwanted 
denominators, including the left-over factor of sinh(gz/) / sinh u from the previous vertex. If 
on the other hand /3 and q have a common factor s then we shall have a left-over factor of 
sinh I// sinh(sz/) besides factors containing Laurent polynomials in y. Furthermore a.i + aj + (3 
will have the same common factor s. The analysis can now be repeated for the next vertex. 
Proceeding this way, and using the fact that the initial charge 7 is taken to be primitive, one 
can prove that at the end all the denominator factors cancel, and we are left with a Laurent 
polynomial in y, proving the desired result. 



C. Equivariant volumes and indices in dipole halos: n = 4, 5 

In this section, we provide explicit results for the equivariant volume and equivariant index 
of the moduli space of n-centered dipole halo configurations with 4 and 5 centers. This serves 
as a check on our minimal modification hypothesis and on the recursion relations derived in 
§6. 3, §6. 4, and provides useful insight of the fixed points responsible for these contributions. 
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C.l Four-centered case with distinct centers 

For = 4 we have to treat several different cases, depending on the value of //2 relative to 
g3, g4 and gs + q^. We label the four possible cases as in [29], see Figure 2. 





Figure 2: Polytopes associated to 4-center dipole halos in ciiarge regimes A, B, C, D respectively, 
in clockwise order starting from top-left corner. 

Case A: + 5'4 < | 

In this case, the contribution to grci{o:i, ■ ■ ■ ,0^4; y) from collinear solutions, given in (6.10), 
is: 

(C.l) 

This has finite y — )■ 1 limit and hence, according to our proposal in §5, should be the complete 
answer. Indeed, it can be checked that (C.l) agrees with the exact result (6.16). 



-1 



iy-y~^Y 
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On the other hand, the equivariant volume evaluates to 

93 



S{I;q3,qi; u) =(-1 







, sinh[(/ - 2(m3 + m4)u] 
dm4 



sinh(/i/) + sinh[(2g3 - /)z/] + sinh[(2g4 - J)z/] + sinh[(/ - 2q3 - 2g4)z/]) 



Replacing z/ — )■ sinh z/ in the denominator we get 



S{I;q3,qi;u) 



3— {sinh(/z/) + sinh[(2g3 - + sinh[(2g4 - 
4 sinh u 



(C.2) 



(C.3) 



+ sinh[(/-2g3-2g4)H} • 



Thus already agrees with the exact result 5* given in (C.l). Hence the prescription below 
(6.17) and the minimal modification hypothesis gives the same result. 

Case B: gg, ^4 < | < gs + 94 

In this case (6.10) leads to the following contribution to (7rcf (o^i, . . . , 04; y)'- 

(-1)^-1 

{y-y 

This diverges in the y — )■ 1 limit. Hence we must add extra contributions of the form given 
in (5.5) to have a Laurent polynomial. This leads to 



y-I ^ ^-7+2,3 + y-I+2,, 



(C.4) 



Sil;q3,q4;y) 



yl _ yl~2,.^ _ yl~2,, _ y^I 



(y-y ^)^ . 

-(2/-r')(2g3 + 2g4-/) 



(y - y-')' 
1 



yl _ yI-2,S _ yI-2,, 



y-I+2qs _^ y-I+2q4, 



for / odd 



y-l y-l + 2q, ^ y-l+2q. 



(C.5) 



-^(2/'-2/"')(2g3 + 2g4-/) 



for / even . 



The need for adding correction terms shows that in this case there are scaling solutions. It 
can be checked that (C.5) indeed agrees with the exact result (6.16). 
On the other hand, the equivariant volume evaluates to 

S{I- gg, g4; = ^ {sinh(/z/) + sinh[(2g3 - + sinh[(2g4 - /)H - (2g3 + 2g4 - W • 

(C.6) 



The last term is recognized as the contribution of the submanifold of fixed points M.^^"^"^ in 
(6.12), with symplectic volume E{I, q^, q^) = (— 1)^^(53 + ^4 — 1/2). Replacing u by sinhz^ in 
the denominator and in the last term of the numerator, we arrive at 



S{r,q3,q4;u) 



(-1) 



7-1 



4 sinh u 



sinh(Ji/) + sinh[(2g3 — J)z/] + sinh[(2g4 — J)z/] 
— (2^3 + 2^4 — /) sinh z/| . 



(C.7) 
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Again, this differs from 5* given in (C.5) by a term that vanishes as z/ — )■ ±00: 



qs, qA\ y) = S{I; q3, qi, y) + (-1 



j-i^qs + 2g4 - J 



: / odd 



4 sinh^ ly \ sinh u — ^ sinh 2z/ : / even j 



(C.8) 

Thus the minimal modification hypothesis and the prescription given below (6.17) agree with 
each other and the exact result. 

Case C: ga < ^ < g4 

In this case (6.10) leads to the following contribution to (7rcf (^i, • • • , ^4; 2/)^ 



Scoii{l; gs, q4; y) 



yl _ _ y-I + ^-/ + 2.3 



(C.9) 



{y-y~^f 

This diverges in the y — )■ 1 limit. Hence we must add extra contributions of the form given 
in (5.5) to have a Laurent polynomial. This leads to 



S{i;q3,qi;y) 



{y - y~'r 



for / odd 



{y-y ^)^ 



/ - y'-''' - y-' + y-'"^^'-' - {y" - y-')q, 

for / even , 



(C.IO) 



in agreement with the exact result (6.16). 

On the other hand, the equivariant volume leads to 



S{r, q3,q4;i^) 



{sinh(/i^) + sinh[(2g3 — I)iy] — 2^31^} 



(C.ll) 



Again, the last term in (C.ll) is recognized as the contribution from the fixed submanifold 
^scai^ with symplectic volume /^'(/, ^3,^4) = (— l)^g3. Following the rules described earlier, 
we get 



S{I;q3,q4;u) = {sm\i{Iu) + sinh[(2g3 - I)iy] -2^3 sinh z/} . 

4smh u 



(C.12) 



It is easy to see that this differs from S given in (C.IO) by terms which vanish as z/ — )■ ±00. 
Case D: ^3,^4 > | 

In this case (6.10) leads to the following contribution to (^ref (oi, • • • , 04; y): 



Scoii{i] q3,qi;y) 



iv-y ^)^ 

This reduces to the case discussed in §5.3. Therefore Eq. (5.9) gives 

r .1)1^1 



(C.13) 



S{i;q3,qA;y) 



(^_^-l)3 



y -y 



-'-Hy-y-') 



y' -y^'-\ iy' - y~') 



for / odd 
for / even . 



(C.14) 



- 54 - 



Again, the results agree with the exact refined index (6.16). 
The equivariant volume gives 



S{I; gg, q^; v) = (sinh(Jz/) - lu) = S'coii(/; gs, Qa] — — , (C.15) 

and hence 

S{I- g3, g4; = ' (sinh(/z/) - / sinh u) . (C.16) 

4sinh u 

Again this differs from the exact result (C.14) by a term that vanishes as — ±00. The last 
term in (C.15) is moreover recognized as the contribution from the fixed submanifold Ai^'^^, 
with symplectic volume E{I ^q^^q^) = (—1)^1/2. We further comment on this contribution 
in §C.3. 

C.2 Five-centered case with distinct centers 

In this section we compute the equivariant volume for 5 distinct centers. Without loss of 
generality, we assume that < ^4 < 55. If in addition 53 + 54 > 55, then we have the ordering 

i) < 53 < 54 < 55 < 53 + 54 < 53 + < 54 + gs < 53 + '?4 + 55 (C.17) 

If instead 55 > 53 + 54, then 

n) < 53 < 54 < 53 + 54 < 55 < 53 + 55 < 54 + 55 < 53 + 54 + 55 (C.18) 

We now split the discussion according to the position of 7/2 relative to these values, starting 
with case i) and then discussing the appropriate change in case ii). The polytopes arising in 
case i) are depicted in Figure 3. 
For 53 + 54 + 55 < //2, 

(-1)^ 

'S'(/; 53, 54, 55; i^) = g^4 (cosh(Ji/) - cosh[(J - 253)//] - cosh[(/ - 254)//] - cosh[(/ - 255)//] 

+ cosh[(/ - 253 - 254)//] + cosh[(/ - 253 - 255)z/] 

+ cosh[(/ — 254 — 255)//] — cosh[(/ — 253 — 254 — 255)1^]) . 



(C.19) 



For 54 + 55 < 1/2 < 53 + 54 + 55, 



(-1)^ 

S{I] 53, 94, 55; ^) = g^4 (cosh(/i^) - cosh[(/ - 253)//] - cosh[(/ - 2qi)u] - cosh[(/ - 255)1^] 
+ cosh[(/ - 253 - 254)//] + cosh[(/ - 253 - 255)//] 

+ cosh[(/ - 254 - 2q,)iy]) - ^-^(/ - 253 - 254 - 255)^ - . 

(C.20) 
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Figure 3: Polytopes associated to 5-center dipole halos with (g3,g4,(/5) = (1,2,4) and varying 
values of /. The vertical, left and right axes correspond to 7715,7713,7714. The dimension associated 
to 771 G with 3 = — Yla''^a is suppressed. 



For ga + gs < < g4 + ^5 



cosh(/z/) — cosh[(/ — 2^3)//] — cosh[(/ — 2^4)//] — cosh[(/ — 2^5)//] 



+ cosh[(/ - 2^3 - 2g4)z/] + cosh[(/ - 2^3 - 2^5)//]] 
+ ^^;f^3(/-g3-2g4-2g5) 



For g3 + '?4 < //2 < Qs + Qb, 



i-iy 



'cosh.{Iu) — cosh[(/ — 2^3)//] — cosh[(/ — 2qi)v\ — cosh[(/ 



cosh[(/ - 2^3 - 2qi)v]] + 



^^\p-Sq^q,-AIq, + Aql) + 



(C.21) 

- 2g5)H 
1)^ 



16i/2 



(C.22) 



For gs < 1/2 < qs + q^, 
S{I]q3,q4,q->]v) 

0/ ■ 

+ 



'cosh(/z/) — cosh[(/ — 2g3)z/] — cosh[(/ — 2q4)v\ — cosh[(/ — 2^5)//]) 



- 2/(g3 + g4 + gs) + 2(g32 + ql + ql)) + 



For g4 < 1/2 < gg, 

5'(/;g3,g4,g5;'^) 



[cosh^Iu) — cosh[(/ — 2g3)z/] — cosh[(/ — 2g4)z/]) 



+ 



^^'(/2-4/(g3 + g4)+4(g32 + g42)) + 



(C.23) 



(C.24) 



16z/2 
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For ga < 1/2 < q^, 

Sil; ga, g4, q,; i^) =^ (cosh(/z/) - cosh[(/ - 2qs)i^]) - ^g3(/ - Qs) • (^-25) 



For < //2 < gg, 



gs, g4, gs; ^y) = ^-^^ [2 cosh(/z/) - u^I^ - 2] . (C.26) 



This case is further discussed in §C.3. 

In case (ii), the region g4 < 1/2 < gs+gs instead sphts in three regions: gs < 1/2 < gs+gs, 
where S is still given by (C.22), ga + g4 < 1/2 < gs, where 

(-1)^ 

S{I; gs, 94, gs; J^) = o. 4 (cosh(Jz/) - cosh[(J - 2q3)u] - cosh[(J - 2qi)iy] 



^ (C.27) 

r 

— n^nA . 

2z/2 



+ cosh[(/ - 2g3 - 2g4)i/]) ^^^q-m 



and g4 < J/2 < ga + g4, where S is still given by (C.24). 
C.3 Multi-equivariant volumes 

To interpret the above results as a sum over isolated and non-isolated fixed points, it is 
useful to compute the equivariant volume for the most general torus action on A^„, 

5(/,{g4;^,K}) = (-l/""^^ L < dma ■ ■ ■ dm. e^^---- ^^^^^^^ ~ ^ ^"=^ ^"^^^ . 

(C.28) 

and compare it to the corresponding equivariant volume of the fixed submanifold M.^^"^"^ 
i5;(/,{ga};{^/a}) = (-l)'-"+' 1 dma---dm„e2S^"™"5( Vm,-^ I . (C.29) 

Jo<ma<ga \ a J 

We shall refer to (C.28) and (C.29) as the multi-equivariant' volume of and Ai^^^, 
respectively. We shall compute these equivariant volumes in two simple cases with n = 4 
and n = 5 centers, which demonstrate that the non-isolated fixed point contribution to 
S{I, {qa}', i^; {^a}) is closely related to E{I, {qa}', i^; {i^a}), though not identical. 
For n = 4 and 1/2 < q3,q4:, Eq. (C.28) evaluates to 



S{I,q3,q4;u,U3,Ui) =(-1 



These four contributions correspond to the four vertices of the polytope V (see Fig. 2, bottom- 
left graph): the first two arise from coUinear configurations (7713 4 = 0, m = ±//2) while the 
last two are of scaling type, with j = (ma = 1/2, = m = or m4 = 1/2, ma = m = 0). 
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Rescaling Ua by a common factor e and taking the limit e — t- (as we shall always do 
when taking the limit Ua — ?► 0), Eq. (C.30) reduces to (C.15). On the other hand, the 
multi-equivariant volume of Ai^^^^ is given by 

E(I,,.,,M = (-1)'-' [^^—^ - . (C.31) 

This differs from the second line in (C.30) for general values of u^, z/4, although it agree with 
it in the limit 1^3,2/4 — > 0, after rescaling by — l/(2z/^). The difference between (C.31) and 
the second line of (C.30) should originate from the Euler class of the normal bundle of A^^"^^^ 
inside Mn, which appears in the denominator of the localization formula. The comparison 
of the two formulae shows that this Euler class should contribute a factor of jiv^ — z/^) at 
each of the fixed points of the toric action. 

To see that such corrections can be important even in the limit — 0, let us consider the 
case n = 5 and //2 < 53, ^4, gs, and identify the fixed points contributing to the equivariant 
volume computed by direct integration in (C.26). In this case, Eq. (C.28) evaluates to 



S{I, qs, 94, gs; i^] 1^3, U5) ={-lY 



e 



+ 



lQiy{iy - iys){iy - i'i){u - z/5) 



16u{u + U3){U + Z/4)(Z/ + Z/g) 8(z/2 - U^^U^ - - Z/5) 



8(z/2 - Z/|)(z/4 - Us){u4, - Z/5) 8(z/2 - Z/|)(z/5 - Z/3)(z/5 - Z/4) 

(C.32) 

These five contributions correspond to the five vertices of the polytope V, displayed on the 
bottom-right corner of Fig. 3 (after restoring the direction along m): the first two arise from 
coUinear configurations (m3^4^5 = 0,m = ±//2) while the last three are of scaling type, with 
j = {rris = //2, 777,4 = ^5 = m = and permutations thereof). Rescaling z/^ by a common 
factor e and taking the limit e — )■ 0, (C.32) reduces to (C.26). In particular, the first two 
terms in (C.32) have a smooth limit at z/^ — and reproduce the first term in (C.26). 
The second, 0{P /v"^) term in (C.26) arises by expanding e*^"^ to second order in z/^, while 
the last, 0{l/v^) term in (C.26) arises by expanding l/(i^^ — z/^) to second order in z/q. In 
contrast, the multi-equivariant volume of A^^*^*^' is given by 



E{l, g3, g4, gs; 1^3, 1^4, ^^s) =(-i)^ 



gl^3^ q^aI qI^sI 

+ 



4(z/3 - Z/4)(z/3 - z/5) 4(z/4 - Z/3)(z/4 - Z/5) 4(z/5 - Z/3) (z/3 - Z/4) 

(C.33) 

This reduces to £^(7,93,^4,95) = (—1)^/^/8 in the limit z/^ 0. Thus, after rescaling by 
a factor — l/(2z/^), the multi-equivariant volume —-£(/, 93, 94, gs; z/3, z/4, z/5)/(2z/^) correctly 
accounts for the 0{P/v'^) term in (C.26), but fails to reproduce the 0(l/z/^). Again, this 
indicates that the Euler class of the normal bundle of A^^'^'*' which appears in the denominator 
of the localization formula should produce an additional factor z/^/(z/^ — z/^) at each fixed 
point. 
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More generally, we expect that the linear combination of equivariant volumes E{I, {qa}) 
appearing in (6.44) can be interpreted as the integral of Ch(£,i^)/ Eu(A^A^^'^'^^) over the 
fixed submanifold Ai^^^^. Similarly, we expect that the analog linear combination of equiv- 
ariant indices E{I, {qa}) which would appear in a similar formula for AS* corresponds to the 
equivariant integral (2.35). It would be interesting to carry this out in detail. 



References 

[1] F. Denef, "Supergravity flows and D-brane stability," JHEP 08 (2000) 050, hep-th/0005049. 

[2] F. Denef, "Quantum quivers and Hafl/hole halos," JHEP 10 (2002) 023, hep-th/0206072. 

[3] B. Bates and F. Denef, "Exact solutions for supersymmetric stationary black hole 
composites," hep-th/0304094. 

[4] J. D. Bekenstein, "Black holes and entropy," Phys. Rev. D7 (1973) 2333-2346. 

[5] S. W. Hawking, "Particle creation by black holes," Commun. Math. Phys. 43 (1975) 199-220. 

[6] S. Ferrara, R. Kallosh, and A. Strominger, "A^ = 2 extremal black holes," Phys. Rev. D52 

(1995) 5412-5416, hep-th/9508072. 

[7] S. Ferrara and R. Kallosh, "Universality of supersymmetric attractors," Phys. Rev. D54 

(1996) 1525-1534, hep-th/9603090. 

[8] S. Ferrara, G. W. Gibbons, and R. Kallosh, "Black holes and critical points in moduli 
space," Nucl. Phys. B500 (1997) 75-93, hep-th/9702103. 

[9] R. M. Wald, "Black hole entropy in the Noether charge," Phys. Rev. D48 (1993) 3427-3431, 
gr-qc/9307038. 

[10] K. Behrndt, G. L. Cardoso, B. de Wit, R. Kallosh, D. Liist, and T. Mohaupt, "Classical and 
quantum n = 2 supersymmetric black holes," Nucl. Phys. B488 (1997) 236-260, 
hep-th/9610105. 

[11] B. de Wit, G. Lopes Cardoso, and T. Mohaupt, "Corrections to macroscopic supersymmetric 
black-hole entropy," Phys. Lett. B451 (1999) 309-316, hep-th/9812082. 

[12] G. L. Cardoso, B. de Wit, J. Kappeli, and T. Mohaupt, "Supersymmetric black hole 
solutions with r**2 interactions," hep-th/0003157. 

[13] A. Sen, "Entropy Function and AdS(2)/CFT(l) Correspondence," JHEP 11 (2008) 075, 
0805.0095. 

[14] R. K. Gupta and A. Sen, "Ads(3)/CFT(2) to Ads(2)/CFT(1)," JHEP 04 (2009) 034, 
0806.0053. 

[15] A. Sen, "Quantum Entropy Function from AdS(2)/CFT(l) Correspondence," Int. J. Mod. 
Phys. A24 (2009) 4225-4244, 0809.3304. 

[16] J. de Boer, S. El-Showk, I. Messamah, and D. Van den Bleeken, "Quantizing N = 2 
Multicenter Solutions," JHEP (2009) 002, 0807.4556. 



- 59 - 



[17] J. Manschot, B. Pioline, and A. Sen, "Wall Crossing from Boltzmann Black Hole Halos," 
1011.1258. 

[18] B. Pioline, "Four ways across the wall," proceedings of the workshop 'Algebra, Geometry 
and Mathematical Physics', Tjarno, Sweden, October 2010, 1103.0261. 

[19] M. Kontsevich and Y. Soibelman, "Stability structures, motivic Donaldson-Thomas 
invariants and cluster transformations," 0811.2435. 

[20] D. Joyce and Y. Song, "A theory of generalized Donaldson-Thomas invariants," 0810.5645. 

[21] J. Manschot, "The Betti numbers of the moduli space of stable sheaves of rank 3 on P^," 
1009.1775. 

[22] J. J. Duistermaat and G. J. Heckman, "On the variation in the cohomology of the symplectic 
form of the reduced phase space," Invent. Math. 69 (1982), no. 2, 259-268. 

[23] F. Denef and G. W. Moore, "Split states, entropy enigmas, holes and halos," 
hep-th/0702146. 

[24] E. Andriyash, F. Denef, D. L. Jafferis, and G. W. Moore, "Wall- crossing from 
supersymmetric galaxies," 1008.0030. 

[25] M. F. Atiyah and R. Bott, "A Lefschetz fixed point formula for elliptic complexes. I," Ann. 
of Math. (2) 86 (1967) 374-407. 

[26] M. F. Atiyah and R. Bott, "A Lefschetz fixed point formula for elliptic complexes. II. 
Applications," Ann. of Math. (2) 88 (1968) 451-491. 

[27] I. Bena, C.-W. Wang, and N. P. Warner, "Mergers and Typical Black Hole Microstates," 
JHEP 11 (2006) 042, hep-th/0608217. 

[28] I. Bena, C.-W. Wang, and N. P. Warner, "Plumbing the Abyss: Black Ring Microstates," 
JHEP 07 (2008) 019, 0706.3786. 

[29] J. de Boer, S. El-Showk, I. Messamah, and D. Van den Bleeken, "A bound on the entropy of 
supergravity?," JHEP 02 (2010) 062, 0906.0011. 

[30] A. Sen, "Arithmetic of Quantum Entropy Function," JHEP 08 (2009) 068, 0903. 1477. 

[31] A. Dabholkar, J. Gomes, S. Murthy, and A. Sen, "Supersymmetric Index from Black Hole 
Entropy," 1009.3226. 

[32] M. Vergne, "Applications of equivariant cohomology," in International Congress of 
Mathematicians. Vol. I, pp. 635-664. Eur. Math. Soc, Ziirich, 2007. 

[33] J. de Boer, F. Denef, S. El-Showk, I. Messamah, and D. Van den Bleeken, "Black hole bound 
states in AdS^ x 52," JHEP 11 (2008) 050, 0802.2257. 

[34] N. Berline, E. Getzler, and M. Vergne, Heat kernels and Dirac operators, vol. 298 of 

Grundlehren der Mathematischen Wissenschaften [Fundamental Principles of Mathematical 
Sciences]. Springer- Verlag, Berlin, 1992. 



- 60 - 



[35] J. J. Duistermaat, The heat kernel Lefschetz fixed point formula for the spin-c Dirac 

operator. Progress in Nonlinear Differential Equations and their Applications, 18. Birkhauser 
Boston Inc., Boston, MA, 1996. 



[36] D. Gaiotto, G. W. Moore, and A. Neitzke, "Framed BPS States," 1006.0146. 

[37] S. Lee and P. Yi, "Framed BPS States, Moduli Dynamics, and Wall-Crossing," 1102.1729. 

[38] N. Berline and M. Vergne, "The equivariant index and Kirillov's character formula," Amer. 
J. Math. 107 (1985), no. 5, 1159-1190. 

[39] S. Alexandrov, D. Persson, and B. Pioline, "Fivebrane instantons, topological wave functions 
and hypermultiplet moduli spaces," 1010.5792. 

[40] M. Shmakova, "Calabi-Yau Black Holes," Phys. Rev. D56 (1997) 540-544, hep-th/9612076. 

[41] B. de Wit, J. Kappeli, G. Lopes Cardoso, and T. Mohaupt, "Examples of stationary BPS 
solutions in A/" = 2 supergravity theories with E? interactions," Fortsch. Phys. 49 (2001) 
557-563, hep-th/0012232. 



- 61 - 



